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Anyons obtained from a finite gauge tlieory have a computational power that depends on the 
symmetry group. The relationship between group structure and computational power is discussed 
in this paper. In particular, it is shown that anyons based on finite groups that are solvable but not 
nilpotent are capable of universal quantum computation. This extends previously published results 
to groups that are smaller, and therefore more practical. Additionally, a new universal gate-set is 
built out of an operation called a probabilistic projection, and a quasi-universal leakage correction 
scheme is discussed. 
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I. INTRODUCTION 

The two main obstacles to building a practical quan- 
tum computer are the decoherence produced by the en- 
vironment and the need for unitary operators of suffi- 
ciently high precision. Topological quantum computa- 
tion provides a way of encoding quantum information 
in non-local observables that are protected from the en- 
vironment, thereby solving the first problem. In some 
instances, the second problem can also be addressed by 
using operations that only depend on topological invari- 
ants. 

Anyons present a concrete realization of the ideas of 
topological quantum computation that may have practi- 
cal implementations. An anyon is a particle that has ex- 
otic quantum statistics, and exists in a two-dimensional 
space. Anyons carry certain topological charges which 
cannot be locally measured or modified, and can there- 
fore be used to store protected quantum information. 
The charges can be detected, though, using two elemen- 
tary operations called braiding and fusion. In the first 
operation, the positions of two anyons in the plane are 
exchanged, causing their world-lines to braid. Because 
clockwise and counterclockwise rotations can be distin- 
guished in two dimensions, braiding can produce results 
more complicated than the usual bosonic and fermionic 
cases. The second operation involves fusing two anyons 
into a single anyon that carries the combined charges 
of the original particles. In both cases, the results only 
depend on the charges and the topological class of the 
paths involved. These operations can be used as a basis 
of error-free gates that manipulate the stored informa- 
tion. 

Many different models of anyons can be constructed 
by specifying different spectra of charges together with 
a set of braiding and fusion rules. A convenient and 
physically inspired set of models can be obtained from 
the electric and magnetic charges of a two-dimensional 
finite-group gauge theory. These models depend on a 
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finite group G, which acts as the symmetry of the gauge 
theory. While every finite group produces a consistent 
model, the computational power of the resulting anyons 
depends on the structure of the group. 

Previous work by the same author 1] has shown that 
finite non-solvable groups produce anyons capable of uni- 
versal quantum computation. However, the smallest fi- 
nite non-solvable group is A^, the even permutations 
of five objects, which has 60 elements. Unfortunately, 
anyons with a large symmetry group are less likely to 
be found in nature, and are also harder to engineer. A 
more desirable symmetry group would be , with only 6 
elements. The purpose of this paper is to study the feasi- 
bility of quantum computation with these smaller groups. 
In fact, it will be shown that the groups that are solv- 
able but not nilpotent, which includes S'3 as the small- 
est case, produce anyons capable of universal quantum 
computation. The caveat, though, is that the construc- 
tions in this paper require both electric and magnetic 
charges, whereas magnetic charges alone were sufficient 
in the non-solvable case Q. The use of electric charges 
complicates the procedure significantly, and will occupy 
the bulk of the discussion. 

The ideas of this paper and its predecessor are built 
on the foundations laid out by Kitaev 01 , who introduced 
the notion of a quantum computer based on anyons. The 
first concrete description for the group A5 was done by 
Ogburn and Preskill in Refs. |^|^. An unpublished con- 
struction for the group S3 was also worked out by Kitaev, 
and its use of electric charges served as a basis for much 
of the present work. 

The organization of this paper is as follows: Section UTI 
contains a review of the basic properties of anyons, and 
describes the notation used in this paper. The next two 
sections prove the universality of anyons based on groups 
that are semidirect products of certain cyclic groups of 
prime order, which includes the important case of S3. 
Section IIIII constructs an abstract set of gates out of 
the fundamental anyon operations, whereas Section IIVI 
proves that this gate-set is universal. In Section the 
discussion is expanded to general finite groups, and the 
relationship between group structure and computational 
power is established. This section will also review the 
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definitions of solvabiHty and nilpotency. The main re- 
sult of this paper, which is the feasibihty of universal 
quantum computation with anyons from groups that are 
solvable but non-nilpotent, is proven in Section IVII The 
discussion in Section IVTI is motivated by Section UTTl and 
includes many of the same steps, but the details are sig- 
nificantly more complicated. Finally, Section IVIII dis- 
cusses a leakage correction scheme that can be applied 
to anyons, as well as many other quantum systems. 



II. REVIEW AND NOTATION 

In this section we present a brief review of the anyon 
properties and operations that will be used throughout 
the paper. Notation for anyon pairs, qudit bases, and a 
special type of gate will be introduced. 

The gauge theory model for anyons that is used in this 
paper was first presented in Ref. and is summarized 
in the review in Ref. |^. Our notation will be closer to 
the one used in the author's previous work though. 
In the interest of space, we present only a very brief re- 
view of the basic ideas and refer the reader to the above 
references for further details. The only new ideas intro- 
duced in this section beyond Ref. are the notation for 
electric charge pairs and a gate called the probabilistic 
projection. 



A. Anyon model and spectrum 

To fully characterize a system with anyons, we must 
specify a set of braiding and fusion rules. A set of con- 
sistent rules can be obtained from the behavior of elec- 
tric and magnetic charges in a finite-group gauge theory. 
Though many other anyon models exist, including models 
with continuous groups, and models with a finite spec- 
trum that cannot be obtained from a gauge theory, only 
the finite-group gauge-theory model will be discussed in 
this paper. 

For each finite group G, there is a gauge theory with 
symmetry group G that contains anyons. The anyonic 
spectrum of the finite-group gauge theory consists of elec- 
tric charges, magnetic charges, and particles called dyons 
which carry both electric and magnetic charge. The mag- 
netic charges, also known as fluxes, are labeled by ele- 
ments of the group G. The electric charges are labeled by 
an irreducible unitary representation R of G, and have 
an internal state that transforms as a vector under R. 
The dyons are labeled by an element of g G C? and a rep- 
resentation of the stabilizer of g. The dyons, however, 
will not play an significant role in this paper. 



B. Magnetic charge pairs 

We begin by discussing the braiding rules for the mag- 
netic charges, which will be the most important particles 



in this paper. The basic rule for magnetic charges is the 
following: when two fiuxes are exchanged, the flux of one 
is conjugated by the flux of the other. Though this is 
the basic interaction that will be used between magnetic 
charges, it has the undesirable consequence that moving 
single fluxes through the system can introduce unwanted 
correlations. Therefore, it will be necessary to work with 
pairs of fluxes of trivial total flux. 

For any g G G, we deflne the state \g) to denote a 
magnetic charge pair, where the first anyon has flux g, 
and the second anyon has flux g^^. Because the effects 
of actions on the compensating flux g^^ will mimic the 
effects on the flux g, we will generally not mention them 
explicitly. In fact, we shall refer to the state \g) as "a 
state of flux g" by which we describe the flux of the first 
anyon, rather than the total flux of the anyon pair which 
shall always be trivial. 

A general state of n magnetic charge pairs has the form 



I*) = 



E 



\9i) 



«> \gn) 



(1) 



where the complex quantum amplitudes. 

Due to the existence of superselection sectors, which will 
be discussed below, the in the above sums are re- 
stricted to a single conjugacy class, which may be dif- 
ferent for each i. 

Because we are dealing with pairs of trivial total flux, 
any two states can be swapped simply by exchanging the 
physical position of the anyons: 



si,S2eG 



9i) ® \92) 



E 



ag^,g^g2)®\9l)■ (2) 



By repeatedly exchanging pairs of adjacent anyons, any 
two pairs of anyons can be exchanged. This operation will 
also allow us to move pre-made ancillas into the compu- 
tational space, and to move anyons that have been fused 
out of the computational space. 

The basic interaction between pairs is a pass-through 
operation by which one pair passes in between a second 
pair. The result of the operation leaves the first pair 
invariant, but conjugates the second pair by either the 
flux of the first pair: 



E ^9ug2\9i) ®\92) 
si,S2eG 



E "si.Sslffl) 
Si,S26G 



1515291 



(3) 



or its inverse: 

E "91,32 l5i> I52) — *E "9i>92l5i> l5r^525i>, 

ffl,S26G 9l,926G 

(4) 

depending on the direction of the pass-though. By using 
the swap operation, the pass-though can be performed 
on any two pairs of anyons. 

Furthermore, the above operation can be generalized 
to a conjugation by a function of the fiuxes of a set of 
anyons. That is, consider a function / : G" G that 
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can be written as a product of its inputs, their inverses, 
and fixed elements of G. For example, 



/ (31,52) ci^i ^C2giC3.giC4.g2C5, 



(5) 



where the {q} are fixed elements of G, and this case has 
n — 2. Then, if we assume the existence of ancillas of 
the form \g) for each g £ G, we can perform the unitary 
transformation 

Iffi) «) I52) «> Iffa) — > 1,91) «> 1.92) «> 1/(51, 52)53/(31, 52)"^), 

(6) 

where we have denoted its action on basis elements, and 
the general transformation follows by linearity. The op- 
eration is implemented by conjugating in sequence by the 
entries of /, starting from right to left, where an ancilla 
of known flux is used for every fixed element of /. In 
general, a conjugation by a function can be performed 
with any state as target, and any number of inputs, as 
long as the function can be written in product form. 



C. Electric charge pairs 



where R{g)M is the matrix obtained by left multiplying 
M by the element g in the representation R. Similarly, 
if we act on the second charge, we obtain 



J®C/(5)|M), 



1 



J2M,,Rl{g)\in) 



\jR') 



= \MR{g-'))^ 



(10) 



where we have used the fact that R is unitary. 

Note that, just as in the case of the magnetic charges, if 
we have a function /({ffi}) of some anyon fluxes, written 
out in product form, then we can apply this function to 
our charges 



\M)^^UU)®I\M)^ = \R{f)M)^ 



(11) 



by applying sequentially from right to left the elements 
of the product. 



In addition to pairs of magnetic charges, this paper will 
often deal with pairs of electric charges, where the first 
charge transforms under the irreducible representation 
i?, and the second charge transforms under the complex 
conjugate representation R* . Of course, for some repre- 
sentations R* ~ i?, which will not be a problem for what 
follows. 

We introduce the bases and {|i_R*)} on which 

the representations act. The indices i, j take values from 
1 to d_R, the dimension of the representation. We assume 
that the basis vectors are compatible in the sense that 

{^RAR*{9)\JR') = {^R\R{g)\3RT■ (7) 

The combined state of the two charges is spanned by 
the vectors \iR)®\jR* ) and can be described by specifying 
a d X d matrix M 

\M)n^^J2MM<^\jR.), (8) 

where we have introduced a convenient normalization fac- 
tor. 

We will be interested in the braiding and fusion prop- 
erties of these states. However, when two electric charges 
move past each other, even when they are not in pairs, 
their charges remain unchanged. It is only the magnetic 
fluxes that have an effect on the electric charges. In par- 
ticular, when a magnetic flux g goes around an electric 
charge, the flux remains invariant, but the charge trans- 
forms as if multiplied by g in the representation R. Start- 
ing with a state |M)jj, if the flux circles the first electric 
charge, then it becomes 

U{g)i»I\M)j^ = -^y^Rik{g)Mkj\iR)®\3R,) 

= \R{9)M)r (9) 



D. Superselection sectors, fusion, and vacuum pairs 

Before describing the fusion rules for the magnetic and 
electric charges, we need to address the issue of supers- 
election sectors, which is familiar to particle physicists. 
A superselection sector is a subspace of a Hilbert space 
that is invariant under all the implementable transforma- 
tions. A useful analogy is to consider the Hilbert space of 
a particle called the nucleon, spanned by the four states 

|OT),|Oi),|iT),|U), (12) 

corresponding to a spin-i particle with two possible 
charge values. This is nothing more than the direct sum 
of the Hilbert spaces of the proton and the neutron: 

'^nucleon ~ ^proton ® ^neutron- (1^) 

At the energies of atomic physics, it is not possible to 
measure in the proton plus neutron basis, or to perform 
a unitary rotation along this direction. Therefore, we 
could say that a nucleon automatically decoheres into 
either a proton or a neutron. 

A similar situation occurs with the anyons. Each 
conjugacy class of G is a magnetic charge superselec- 
tion sector. The irreducible representations are the elec- 
tric charge superselection sectors. When given an un- 
known anyon — for example, an anyon created from the 
vacuum — we can assume that it has decohered into a spe- 
cific, though possibly unknown, conjugacy class and/or 
irreducible representation. Furthermore, when storing 
quantum information, it will be important to keep the 
computational space in a single superselection sector to 
avoid decoherence. 

Let I'l') be a pair of anyons created from the vacuum. 
We may assume that each anyon has decohered into a 
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specific superselection sector. Furthermore, because a 
vacuum pair must consist of a particle with its antipar- 
ticle, the two superselection sectors are related. That 
is, the pair must have vacuum quantum numbers and be 
able to fuse back into the vacuum. Therefore, if the first 
anyon is a magnetic charge with flux in a given conjugacy 
class, the second anyon will be a magnetic charge with 
flux in the inverse conjugacy class. If the first anyon is 
an electric charge of representation R, then the second 
anyon will be a electric charge of the complex conjugate 
representation. Finally, if one anyon is a dyon, then so is 
the other. 

In the case of magnetic charges, there is exactly one 
state with vacuum quantum numbers in each conjugacy 
class. The state is 



|Vac(C)) 



1 



(14) 



gee 



where C is a conjugacy class of G. Note that, given our 
notation, the above state is an entangled state of two 
anyons. In the case of electric charges, the vacuum state 
for representation R is simply |i?(/))^, where R{I) is the 
dfi, X dR identity matrix. 

The operation of fusion is in a sense the inverse of vac- 
uum pair creation. Fusing two anyons produces a single 
anyon that must carry the total magnetic and electric 
charges of the pair. In the special case when both total 
charges are trivial (i.e., one of the above vacuum states) 
the state can fuse into the vacuum, leaving no particle be- 
hind, and transferring its energy to some other medium 
such as photons. In theory, this case can easily be de- 
tected in the laboratory, and is the primary way of ob- 
taining measurement results. 

In the case of magnetic charges, the net resulting flux 
is just the product of the two fluxes, where the order- 
ing of the product depends on some conventions which 
will not be important here. While one of our standard 
anyon pairs always has trivial total flux, we sometimes 
may fuse anyons from different pairs to determine if their 
flux is equal. Even if the total flux is trivial, though, 
the pair may not fuse into the vacuum but may pro- 
duce an electric charge. This will be the case if the state 
transforms non-trivially under simultaneous conjugation 
of both anyons. 

The fusion of two electric charges can only produce an- 
other electric charge (or the vacuum, which is the charge 
carrying the trivial representation). To calculate the pos- 
sible products of fusion, note that fusion implies that a 
flux can no longer be braided around only one of the two 
electric charges. Mathematically, it is a restriction to the 
diagonal transformations 



\M)j 



U{g) ® U{g)\M)^ = \R{g)MR{g-')) (15) 



However, the above action of the group is not irre- 
ducible on this space. The vector space spanned by 
all possible states |M)^ decomposes into invariant sub- 
spaces. The invariant subspaces correspond to electric 



charges transforming under irreducible representations. 
The probability of obtaining each irreducible representa- 
tion corresponds to the magnitude of the state vector pro- 
jected down to the appropriate invariant subspace. Fur- 
thermore, after fusion, it is no longer possible to measure 
the relative phase between the different representations 
and therefore decoherence occurs in the representation 
basis. 

The net result of fusion is a mixed state of different rep- 
resentations. Which representations occur is determined 
by the decomposition of R{g)iS)R* (g) into irreducible rep- 
resentations. The probability of obtaining each of these 
representations is determined by the projection of M to 
the different invariant subspaces. 

In particular, the trace of M is the unique invariant 
under conjugation by G (which is the content of Schur's 
lemma). Therefore the probability of fusion into the vac- 
uum is 



vac 



|(i?(l)|M), 



Tr (M) 



(16) 



E. Requirements for the physical system 

To complete our review of the properties of anyons, 
we will list the operations, ancillas, and measurements 
that we assume are available on any realistic system, and 
which we will use to build our quantum gate-set: 

1. We can braid or exchange any two particles. 

2. We can fuse a pair of anyons and detect whether 
there is a particle left behind or whether they had 
vacuum quantum numbers. 

3. We can produce a pair of anyons in a state that is 
chosen at random from the two particle subspace 
that has vacuum quantum numbers. 

4. We have a supply of ancillas of the form \g) for any 
geG. 

5. We have a supply of ancillas of the form |i?(/))^ 
for any irreducible unitary representation R. 

The last two requirements are the only questionable 
ones, as it is not obvious how to produce this reservoir of 
calibrated electric and magnetic charges. In fact, since 
many of these ancillas will be destroyed during fusion, 
the reservoir will have to have a large number of ancillas 
of each type. 

One of the main difference between the constructions 
in this paper, and the one used in producing computa- 
tions with non-solvable groups 1], is that the latter case 
rec^uired no electric charge ancillas, which may be harder 
to produce. Additionally, Ref. presented a protocol 
for producing the magnetic ancillas for a simple non- 
abelian group. The production of calibrated flux and 
charge ancillas for the groups discussed in the present 
paper, though similar, will not be addressed here. 
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A final note is that the requirement of cahbrated mag- 
netic charge ancillas wiU have to be shghtly modified in 
Section IVl CI in order to work with certain large groups. 



F. Notation for Qudits 

Throughout this paper it will be useful to perform 
computations with qudits rather than the usual qubits. 
We define our computational basis as the states \i) for 
< i < d, where we will assume that d is prime. The 
unitary Z and X gates can be defined as follows: 



XH) = U + 



(17) 
(18) 



where w is a fixed non-trivial d*^ root of unity, and sums 
are understood to be modulo d. As usual, the eigenstates 
of Z correspond to the computational basis. We can also 
introduce the eigenstates of X: 



Td^' 



(19) 



which have the following transformations under the ac- 
tion of our unitary gates: 



Z\i) = 
X\i) = 



(20) 
(21) 



Note that when appropriate, we shall assume all oper- 
ations are modulo d without further comment. 



Furthermore, we demand that if outcome is obtained 
when applying the operation to a state \^), we obtain 
the state 



l*o> 



(23) 



On the other hand, if we get the result 1, we will consider 
the state damaged, and trace it out of our computational 
system. 

As an example consider 



1 



(24) 



and let K, = {|0)}. Applying a probabilistic projection 
to the first qubit, we obtain with probability pppjl the 
state 

l*o) = |0)® |1), (25) 
and with probability l—ppp/2 we obtain the mixed state 



Pi 



PPP 



(1-PPP)|1)(1| + 



(26) 



where we have already traced out the first qubit. Notice 
that if the probabilistic projection onto |0) is applied to 
both qubits simultaneously, it is possible to obtain the 
result 1 twice, but it is not possible to obtain the result 
twice. 



III. BASE CASE: G = xe 



G. Probabilistic projection onto K, 

To conclude with the introduction of notation, we de- 
fine a new type of gate called a probabilistic projection 
onto a subspace. The operation is essentially a projec- 
tive measurement that distinguishes between a subspace 
/C, and its orthogonal complement. However, the opera- 
tion has a one-sided probability of error, corresponding 
to a failure to notice the projection into K.. 

For example, consider an operation that emits a pho- 
ton if and only if the state is projected into the subspace 
JC. The photon is then received at a photodetector that 
has a probability < p < 1 of absorbing the photon. A 
photon will never be detected if the state was projected 
into the complement of /C, but even if the measurement 
projected into /C, the photodetector may remain silent. 

To formalize the idea of a probabilistic projection, let 
/C be a subspace of a Hilbert space Ti, and let P/c be the 
projection onto JC. We define a probabilistic projection 
onto /C as a two-outcome POVM with operators 

Fo = PPPP^, Fi = l- pppPk, (22) 

where < ppp < 1. We say that we can do a probabilis- 
tic projection onto /C if we can do the above operation 
for any fixed ppp. 



Before tackling the general case of groups that are solv- 
able but not nilpotent, we will describe the procedure for 
producing quantum computation using a special type of 
group based on the semidirect product. The construc- 
tion for these groups is very similar to the general case, 
but can be described in more concrete terms. In particu- 
lar, these groups are very useful in eliminating operations 
whose usefulness is unclear in the general case, but that 
have no computational power when reduced to this spe- 
cial case. 



A. Algebraic structure 

We will be interested in the groups G = Zp xgZq, the 
semidirect product of the cyclic groups of order p and 
q. We assume that p ^ q are both prime and that the 
function Q is non-trivial, which guarantees that G is not 
nilpotent. 

The group can be described using two generators a and 
h which satisfy the relations: 

af^l, 6« = 1, hah-^=a\ (27) 

where specifying an integer t between and p is equiv- 
alent to specifying the function Q : Iq ^ Aut(Zp) used 
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for the semidirect product. We will require that t ^ 1 
which is equivalent to 9 being non-trivial. Furthermore, 
consistency requires that 



= 1 mod p, 



(28) 



which can always be solved for some t as long as q divides 
p — 1. We henceforth assume that p, q, and t have been 
chosen in a self-consistent fashion. 

The best example of one of these groups, and in fact 
the smallest non-abelian group, is S3. This group can be 
expressed as Z3 XeZ2, with t = 2. We can choose a to be 
any order three element such as (123), and we can choose 
6 to be any order two element such as (12). 

The first example of such a group with odd order is 
Z7 X 0Z3 with t = 2 or t = 4, both of which are equivalent. 
One of the most important features of this example is 
that not all the non-trivial powers of a are conjugate 
to one another. The elements a, a}, and form one 
conjugacy class, whereas the elements a'^, a^, and form 
another. 

Both of the above examples will be revisited when 
we discuss group representations and fusion of electric 
charges. 



B. Computational basis 



We choose a qudit computational basis 



(29) 



for < i < p. Note that all these states are unique 
because a^ba~^ = a^'-^~*'>b, and o^~* is a non-trivial gen- 
erator of the group Zp. We are therefore using a complete 
conjugacy class for the computational subspace. 

While the above choice of computational subspace may 
seem arbitrary, most other choices are either equivalent 
or less powerful. The conjugacy classes a^Va~'', for dif- 
ferent non-trivial values of j, are all equivalent. Dyons 
with these fluxes are also equivalent since they are just 
the combination of the above states with electric charges 
that cannot be detected by braiding. Finally, the powers 
of a and pure electric charges are suboptimal as they are 
difficult to entangle (for more on this see the discussion 
on using nilpotent groups in Section lv|) . 

Initializing a quantum computer in this basis is easy, 
as we have assumed the existence of flux ancillas in the 
state |0), which can be used as computational anyons. 
We therefore turn to the task of implementing gates on 
this space. 



C. Operations involving braiding fluxes 

We begin by characterizing the operations that can be 
achieved by braiding fluxes. Fix a target qudit which we 
will be conjugating, and assume that it is in the compu- 
tational subspace. We can conjugate this qudit by the 



fluxes of arbitrary ancillas in the group. It can also be 
conjugated by the fluxes of other qudits, which we will 
also assume to have a definite flux in the computational 
subspace (as the effect of a superposition of fluxes can be 
inferred by linearity). 

Let us begin with the case when only one qudit (in 
addition to the target) is involved. If the source qudit is 
in a state \g), then the target will get conjugated by an 
expression 



fig) = Cigc2gc3 . 



(30) 



for some n, where the {cj} are fixed elements of G corre- 
sponding to the ancillas used. Of course, these elements 
represent the product of any ancillas that were used in 
series, and can equal the identity if no ancillas were used. 

Because of the structure of the group, all the fixed el- 
ements can be expanded as Ci = a^^b^^ for some integers 
ji, ki. Furthermore, since the source flux is in the com- 
putational basis, it can be written out as g = a^6a^^ — 
a for some X. Inserting these expressions, we get 



fig) 



-^^ba^^b''^ ...a^-b''". 



(31) 



Using the group relation 6a* 



*6, we can move all 



the 6's to the right, and combine factors to get 



fig) = a"a'3-6^ 



(32) 



for some integers a, (3 and 5. The effect of each of these 
factors can be considered separately. Conjugating by a" 
is just the application of the gate X". Conjugating by 
a'^^ is just a controUed-X from the source to the target, 
repeated (3 times. Finally, conjugating by b maps \i) to 
I it). This operation can be generated using a controUed- 
X gate, and an ancilla |0): 



|0) 



X- 



|0) 



it) 



where —\/t is computed modulo p. Following the above 
circuit, we can either replace the original qudit with the 
ancilla, or use a swap, which can also be built out of 
controUed-X gates. 

So far we have shown that the X and controUed-X 
gates generate the set of operations achieved by conju- 
gations. However, we have yet to show that these opera- 
tions are in fact included in the set of achievable opera- 
tions. The X gate is rather trivial as it is a conjugation 
by an ancilla of flux a. The controUed-X is a conjugation 
by the function 



fig) = igh 

= (a-^^-%b-^) 



]^-^ 1/(1— ■/:) mod p 

l/{l-t) mod p 



(33) 
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where 1/(1 — t) can be computed modulo p because we 
assumed 1 <t < p. 

The case involving many source qudits, all of which 
can be used to conjugate the target, is very similar to 
the above. The expression can be simplified by moving 
all the 6's to the left, and combining similar factors. In 
the end, the net effect will again be a series of X and 
controlled-X gates. 

Finally, one may wonder about using an ancilla as an 
intermediate step. That is, first we take an ancilla (say, 
g'), conjugate it by some function (say, /) of some qudits, 
and then conjugate the target by the ancilla. However, 
the same effect can be achieved by conjugating the target 
first by f~^, then by g' and finally by /. This procedure 
therefore provides no extra computational power. 

The conclusion is that the operations achievable from 
braiding magnetic charges are exactly those generated by 
the X and controUed-X gates. In fact, the X gate is re- 
dundant as we have assumed the existence of 1 1) ancillas, 
which can be used as control qudits in a controlled-X. 



D. Operations involving fusion of fluxes 

Now we turn to the operations achieved by the fusion 
of magnetic fluxes. For these operations it will be suffi- 
cient to determine whether the two particles fused into 
the vacuum or not, thereby obtaining at most one bit of 
information from each fusion. 

At this point we remind the reader that standard states 
consist of pairs of anyons, whose total flux is trivial. That 
is, the state \g) describes an anyon of flux g paired with an 
anyon of flux g~^. There are therefore two basic choices 
for fusion: we can fuse the two anyons that compose a 
single pair with each other, or we can fuse one of them 
with an anyon from another pair, typically an ancilla. To 
avoid confusion, in the latter case we will always use the 
anyon of flux g (rather than g~^) for the fusion. 

The case of fusion with an ancilla will lead to a mea- 
surement in the Z basis. The fusion of anyons from the 
same pair will lead to a measurement in the X basis. 
However, we will delay the construction of the actual 
measurement gates until the next section. For this sec- 
tion, we will simply describe the fusions as abstract op- 
erations on the computational space by employing the 
construction of probabilistic projections. 

The fusion of an anyon from a state l^*) with an anyon 
ancilla of flux b^^ is a probabilistic projection onto the 
subspace /C = {|0)}. That is, an anyon of flux a^ba~'^ 
can only fuse into the vacuum with a flux if i = 
(modulo p as usual) . When i > there must be an anyon 
left over to carry the non-trivial total flux. When i = 
the fusion can either produce the vacuum state or an 
anyon with non-trivial charge. The probability for fusion 
into the vacuum in this case is 1/p. Furthermore, if we 
fuse into the vacuum we can replace the state with a |0) 
ancilla. Therefore the whole operation is a probabilistic 
projection onto |0) with ppp = 1/p. 



The fusion of two anyons from the same pair is a proba- 
bilistic projection onto the subspace /C = {|0)}. Because 
the total magnetic flux of the pair is always trivial, the 
fusion product must be an electric charge. The charge 
corresponds to a representation of G given by the action 
of conjugation on the anyon fluxes. The state |0) trans- 
forms trivially and corresponds to the vacuum, whereas 
the states for i > 0, are orthogonal to the vacuum and 
correspond to non-trivial representations. In fact, this 
procedure is a probabilistic projection with ppp = 1. 
However, since the state is destroyed during fusion, to 
complete the projection we must be able to produce |0) 
states. This will be discussed below. 

The other choices for fusion are equivalent to a com- 
bination of one of the above measurements and an X 
or controlled-X gate. Fusing with a flux of the form 
a*6~^a~* is equivalent to flrst applying a gate and 
then performing a fusion with b~^. A fusion with any 
other flux can never produce the vacuum if the qudit is 
in the computational subspace. Finally, one can consider 
fusion of anyons from two different qudits. If the state 
of the two qudits is |z) (g) the fusion will only produce 
the vacuum state if i = j. Therefore, the operation can 
be simulated by a controUed-X^^, followed by the fusion 
of the target with a b^^ flux. 

The conclusion so far is that fusion of magnetic charges 
provides us with two new operations: the probabilistic 
projections onto the subspaces |0) and |0), which will 
eventually become measurements in the Z and X bases. 
The only operation that has not been considered is using 
the products of fusion for further operations or fusions. 
This subject will be briefly touched upon after discussing 
fusion of electric charges. 

1. Production of |0) states 

To conclude the discussion on fusion of fluxes, we 
present the construction of |0) states, which were needed 
to complete the probabilistic projection onto |0). 

Just as the state |0) naturally fuses into the vacuum, 
it is also naturally produced from the vacuum. Unfor- 
tunately, producing a pair of anyons from the vacuum is 
just as likely to produce the vacuum state for one of the 
other superselection sectors as it is to produce the state 
|0). Therefore, after producing a vacuum state we must 
measure its superselection sector. Vacuum pairs that are 
produced in the computation subspace (magnetic charge 
in the conjugacy class of 6) will be kept as |0) states, and 
the rest will be tossed out. 

Since measurements are done by fusion, which is a de- 
structive procedure, we must copy the vacuum state be- 
fore measuring the conjugacy class. The procedure starts 
with a pair created from the vacuum and a |0) ancilla: 

I Vac) ® |0), (34) 

and applies to it a swap, made out of the conjugation- 
based controUed-X: 
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they can be easily described in their natural basis as 
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where the circuit depicts the resuh for the case when the 
vacuum pair was created in the computational superse- 
lection sector, in which case |Vac) = |0). 

In the case when the vacuum state was not created in 
the computational supersclcction sector, then the effect 
of the conjugations will be different. However, since the 
conjugations are performed using braiding, which never 
changes the superselection sector, the vacuum state can 
only be transformed into a state that is orthogonal to 
|0) = 

After applying the above controUed-X gates, wc at- 
tempt to fuse an anyon from what was the vacuum state 
with an ancilla of flux b~^. If they fuse into the vac- 
uum, this implies that the vacuum state was created in 
the computational superselection sector, and the above 
circuit worked correctly. The ancillas |0) will have been 
transformed properly into a |0) ancilla, which can be used 
for computation. In the case when the fusion does not 
produce a vacuum state, the swap probably did not pro- 
duced the desired state, so we discard it and start over. 

To summarize, we now have a source of |0) ancillas, 
which can be used as the last step needed to complete 
the probabilistic projection onto |0). 



E. Representations and fusion of electric charges 

Thus far, we have only considered operations involving 
magnetic fluxes. These operations led to a controUed-X, 
and measurements in the X and Z bases. However, these 
gates do not form a universal gate-set. Wc must therefore 
consider operations involving electric charges as well. 

The electric charges transform as irreducible represen- 
tation of the group G. To obtain the spectrum of electric 
charges, as well as their braiding and fusion rules, we 
must therefore discuss the representation theory of G. 

It is easy to sec that the commutator subgroup G" of 
groups of the form G = Zp xg Zq is just G' = Zp. The 
representation theory of G can be obtained by inducing 
representations from G". Starting from the trivial rep- 
resentation on G' , the induced representations are the 
one-dimensional representations where a — > 1 and 6 is a 
g*^ root of unity. 

The rest of the irreducible representations have dimen- 
sion q and are obtained by inducing from the non-trivial 

representations of Zp. The induced representations are 
all irreducible though not necessarily distinct. In fact. 



where oj is the p*'' root of unity of the representation from 
which we are inducing. The matrix for a is diagonal, 
whereas b is the permutation matrix with entries 1 above 
the diagonal. 

Even though the representation theory for these par- 
ticular groups is easy, we will use abstract language to 
describe the fusion rules, which will make the connection 
to the general case clearer. 

Take any non-abelian irreducible representation, and 
consider a pair of electric charges in the state |i?(a*))^. 
What representations do we get if we fuse the two 
charges? The product of fusion is invariant under the 
action of a: 



U{a)(S)U{a)\R{a'))^ 



\Ria)R{a')R{a-^))j^ 
\Rial)n (35) 



and therefore represents the commutator subgroup G' 
by the identity. This implies that the representation is 
abelian! In particular, it is easy to see that the one- 
dimensional subspaces 



|[7^])^= |diag(7^7'^•••,7«^■))^ 



(36) 



with 7"? = 1 are the spaces corresponding to the repre- 
sentations a ^ 1, & ^ 7-^ . 

We will be interested in the quantum amplitude that a 
state |i?(a*))jj fuses into the b 7-^ representation. This 
quantity will be denoted by the fusion amplitude 
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(k-1) 



(37) 



k=l 



with < i < p and < j < q. 

Let l^*) be an arbitrary state entangled with an electric 
charge pair 



i=0 



(38) 



where the |\E'i) denotes (unnormahzed) states of the rest 
of the system. The fusion amplitudes allow \^) to be 
rewritten as 



i=0 j=0 



(39) 



The basis | [7-'])^ labels the total charge of the two anyons 

that comprise the electric charge pair. A fusion of the 
two electric charges, followed by a measurement of the 
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resulting fusion product, will be a measurement in this 
basis. 

Note that the basis |[7-'])7j only spans the diagonal 
subspace of However, this is the subspace con- 

taining all the states |i?(a'))^. The subspaces spanned 
by |i?(6-'a'))^, for some fixed j > 0, are mapped un- 
changed into the space of a single higher-dimensional ir- 
reducible representation, and are therefore not useful for 
the purposes of this paper. 

While the representation R does not appear explicitly 
in the fusion coefficients, it enters implicitly in the above 
expression as the choice for p*^ root of unity uj. Though 
we could use the notation uju, this will not be necessary as 
we will generally work with only one higher-dimensional 
irreducible representation. 

The most important feature of the Fi^j coefficients is 
that \R{a^))j^ is the vacuum state and therefore 



(40) 



which can be verified by direct calculation. Another im- 
portant property is that 



(41) 



for all i > 0. The proof involves showing that a linear 
relation of roots of unity only vanishes if it is a combi- 
nation of the obvious regular polygon relations (which is 
proven in 0). 

A final interesting property is that 



(42) 



which is a consequence of 



|i?(a"'=))^ = \R{b'a^b-'))^ 

= U{b'')(E>U{b'')\R{a'))jj. (43) 

F. Examples 

1. Ss 

The group S3 has three irreducible representations, the 
trivial (identity) representation (where a ^ 1, 6 ^ 1), 
the sign of the permutation (where a — > 1, 6 — s- —1) and 
a two-dimensional one: 



CJ 
U! 



V 

1 01 ' 



(44) 



where w is a non-trivial cube root of unity. The fusion 
amplitudes are 



^0^0 — 1, Fi^o 



1 

2' 



F. 



Fo^i = 0, Fi^i = -z^, F2^i^i^. (45) 

The best way to visualize these coefficients is to 
start with a state |0), and a pair of electric charges in 



the vacuum state of the two-dimensional representation: 
|i?(/))^. Then entangle with a controUed-sum to get 



^Y.\^)\R{a^)), 



-y 



\J) 







(46) 



Fusion of the electric charge pair produces either the 
vacuum (trivial representation) or a charge transforming 
under the sign representation. The probability of getting 
each is 



p — 

^ vac 



^ sgn 



E 
E 





2 


1 






2' 




2 


1 






2' 



(47) 



and the state of the magnetic charges afterwards is one 
of 



1 



I* 



sgn/ 



^^(2|0)-|l)- 



|2)). 



(48) 



These are obtained by multiplying the initial state by 
the appropriate F coefficients and renormalizing to unit 
magnitude. In the case of the second state we also intro- 
duced an extra global phase of i, which is related to the 
arbitrary choice of phase of the |[7-'])jj. states. 



2. 



The group ZiXglj^ has five irreducible representations. 
Three of them are one dimensional and set a — > 1 and b to 
a cube root of unity. The other two are three dimensional 
and are complex conjugates of each other. 

The main new feature of this group is that the non- 
trivial powers of a are not all conjugate to one another. 
This leads to more complicated fusion coefficients. For 
example: 



Fo- 



0, 



Fi^i = |a, Fs^i = ^B, Fe^i = is, (49) 



with 

A = 7^CJ + ^uj'^ + Lo'^ 

B = 7^cj"^ +7cj"^ 

= e^^^sT +e2"5T +6^"*, (50) 

where we have chosen 7 = e^'^^/'^ and lj = ^'^^1'^ . Notice 
how A is close in magnitude to 3 whereas B is close in 
magnitude to 1. 
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G. Operations involving electric charges 

Now it is time to apply the discussion in the previ- 
ous subsections to build a useful operation out of electric 
charges. While there seems to be a wealth of strange an- 
cillas that could be produced using electric charges, most 
of them have complicated relative amplitudes or phases 
that are hard to use in a constructive proof of univer- 
sal computation. We will therefore focus our attention 
on producing an operation that arises naturally from the 
fusion amplitudes: the projection onto the subspace or- 
thogonal to |0). 

Consider a qudit in the state 

l*>=E'/'d^>, (51) 

i=Q 

where the coefficients {ipi} could either be complex num- 
bers, or could represent the state of the rest of the system 
if the qudit is entangled with other qudits. 

We append to the qudit an electric charge pair \R{I))j^ 
in the vacuum state of a non-abelian representation R. 
Using braiding, we can right multiply the state of the 
electric charge by some function / of the qudits flux: 

14-) ® \R{I))j, ^ ® (52) 

i=0 

We have shown in Section Fill CI that the most general 
function is of the form f(i) — a^a^'^b^ . Choosing S ^ 
turns out not to be useful, and choosing a 7^ can be 
used to get projections to the spaces orthogonal to \i) for 
i > 0, but this can be achieved as well with an X gate. 
We will therefore focus on f{i) = a^"^ so that we obtain 
the state 

p— 1 p^i 

E ® \Ria^'))j, -EE PP^-^M ® I b'])R- 
1=0 1=0 i=o 

(53) 

A fusion of the electric charge pair, followed by a mea- 
surement of the resulting electric charge (the feasibility 
of which will be the subject of Section llll HI below), leads 
to a state that is proportional to 

p-i 

J2Fp^^J4'^\^), (54) 
i=0 

where j now labels the result of the measurement in the 
basis I [7^])^. 

Because of the property Fq^j = Sj^, if the measure- 
ment result is j 7^ 0, we will have projected into the space 
orthogonal to |0). Unfortunately, we will have also intro- 
duced undesired relative phases and amplitudes. The 
trick will be to balance these out. 

Consider repeating the above procedure p — 1 times, 
with P taking values from 1 to p— 1. Furthermore, assume 



that in each case the fusion results in j = 1 . The resulting 
state will be, up to normalization: 

p-i /p-i \ p-i 

«E n^/3*-ip''i*)«E'/''i^)' (55) 

i=i y/3=i J i=i 

where we have used the fact that multiplication by i, 
modulo p, is just a rearrangement of the values of (3. 

The above procedure is a probabilistic projection onto 
/C = |0>-^. As usual, if we do not obtain j = 1 as the 
result of each measurement, we just discard the state 
being projected. 

What is the probability of success of the above proce- 
dure? The probability for obtaining j = 1 on the first 
try is 

p— 1 p^i 
P,=i = E l^-iV^d' > mm (iF.^il') E l^^l' ■ (56) 

i=l 4=1 

On subsequent measurements, the state has previously 
been projected to |0)^ and renormalized. Therefore the 
probability of success for each trial is simply bounded by 

P,=i>mm(|i^,^i|'). (57) 

The total probability of success is just the product of 
these quantities. In particular, the probability ppp asso- 
ciated with the probabilistic projection can be bounded 
by 

PPP >min(|F,^i I')'' ' >0, (58) 
j>0 \ / 

where we used the fact that > for i > 0. 

Of course, the above is a underestimation of the prob- 
ability of obtaining a good projection. For example, if 
all the results j were equal to some fixed j > 1 , the same 
argument would show that a correct projection was ob- 
tained. Furthermore, there are many other ways in which 
the relative phases and amplitudes can cancel out. A 
classical computer, with knowledge of the values of Fi^j , 
can keep repeating the procedure until such a cancella- 
tion occurs. The computer would also be required to stop 
after a long sequence of j = results, in which case the 
state would have been projected onto |0). 

In the end, as long as ppp is fixed and finite, we have 
produced the desired probabilistic projection to the space 
|0)^. Different values of ppp will just affect the complex- 
ity of an algorithm as a multiplicative constant. Further- 
more, for the small groups that are likely to appear in the 
laboratory, ppp should be reasonably large. For example, 
in the case of G = 6*3, ppp can be made exponentially 
close to one in the number of measurements. 

It should be noted that because we are working with 
qudits of dimension d = p, and the semidirect product 
requires p > q > 2, the above projection will always be a 
non-trivial operation. In fact, it will always be powerful 
enough to complete a universal gate-set. 
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At this point, all that remains to be done is to prove 
the universality of the gates constructed from the ba- 
sic anyon operations. This will be the subject of Sec- 
tion IIVI However, before closing this section, we shall 
discuss some issues regarding the measurability of elec- 
tric charges, and look at some alternative operations that 
could have been employed. 



H. On the measurement of one-dimensional 
representations 

The feasibility and accuracy of the probabilistic pro- 
jection onto |0)''~ depend crucially on being able to iden- 
tify electric charges carrying one-dimensional represen- 
tations. However, these charges have a special property 
that makes them hard to identify: when only using braid- 
ing, a one-dimensional representation is indistinguishable 
from the vacuum! 

The reason behind the above difficulty is that one- 
dimensional representations of a group G are constant 
on conjugacy classes of G. Therefore, a magnetic charge 
that is braided around one of these electric charges will 
have its state change by an overall phase. These global 
phases are not measurable in quantum mechanics. 

Of course, an interference experiment would produce a 
measurement of the charge. The standard double-slit ex- 
periment, with the electric charge located in between the 
slits, will produce a pattern on the screen that depends 
on the representation of the electric charge. However, 
during the experiment, the anyon will be in a superposi- 
tion of spatial positions which is no longer protected from 
decoherence by topology. Since the interference experi- 
ment can be repeated many times without affecting the 
electric charge, this may not necessarily be a problem. 
However, it does involve working in a regime where the 
anyons can be treated as waves rather than particles. 

On the bright side, these electric charges can also 
be detected by fusion, assuming the availability of elec- 
tric charge ancillas with one-dimensional representations. 
Their fusion rules are particularly simple because these 
states have a one-dimensional internal Hilbert space. 
Furthermore, their fusions always produce unique results. 
If 7(3) and ^' {g) are two one-dimensional representations 
of a group G, then the fusion of the electric charges car- 
rying these representations produces a charge of repre- 
sentation 7"((?) = 7(3)7' (3). A charge will only fuse into 
the vacuum when fused with its conjugate representation. 
Therefore, after a series of fusions that end up producing 
the vacuum state, we can determine the representation 
of the original electric charge. 

In fact, for groups with q = 2 such as 6*3, there is 
a further simplification. In these groups there are only 
two one-dimensional representations: the vacuum and 
the sign representation. Since the fusion of |i?(a*))^ pro- 
duces a one-dimensional charge, if it does not fuse into 
the vacuum, then it must have produced the sign charge. 
Therefore, for these groups, we do not even require one- 



dimensional electric charge ancillas. 



I. Other possibilities 

In this section, we will briefly discuss one last possibil- 
ity for producing useful operations: using the products 
of fusion. Though not strictly needed to complete a uni- 
versal gate-set, this subsection is an interesting study of 
alternative operations and the effects of decoherence dur- 
ing fusion. 

At flrst sight, it appears that the projection onto |0) 
can be done without using electric charges with the fol- 
lowing procedure: first fuse one anyon from the state to 
be measured with a flux. Only the |0) state can fuse 
into the vacuum. If an anyon remains, fuse again with a 
h flux to restore it to its previous state, and pair it with 
its old partner. Repeating the procedure multiple times 
(because the |0) could turn into an electric charge rather 
than the vacuum) yields the desired projection. 

There are, however, two problems with the above con- 
struction. The flrst, and smaller, problem is that when 
fusing with b^^ or with h we could be turning our mag- 
netic charges into dyons. For groups of the form Zp XgZq 
the dyonic electric charges are all one dimensional, how- 
ever, and will therefore have no effect on braiding, as 
discussed in the previous section. The probabilities of 
fusion into the vacuum will be reduced, and therefore so 
will the respective projection probabilities, but they will 
still remain non-zero. In fact, a careful examination of 
the operations constructed so far shows that they work 
with a probabilistic mixture of dyons and regular mag- 
netic charges. 

The second and larger problem, though, is decoher- 
ence. The fluxes a^ba~^b~^ = a**^^~*^ belong, in general, 
to different conjugacy classes and therefore different su- 
perselection sectors. When the quantum state is encoded 
in this form, it is susceptible to decohere into the different 
superselection sectors. 

When does this decoherence occur? It occurs during 
fusion. In general, fusion takes two n-dimensional Hilbert 
spaces H and maps them to one: Hi x 7^2 H3. But 
quantum mechanics is unitary; therefore, what must re- 
ally be happening is a mapping to a tensor product of Ti. 
and the environment: Hi x H2 H3 x £. When two 
states are mapped onto new states that are orthogonal 
in the environment subspace, decoherence occurs. 

How do we know if states will have orthogonal envi- 
ronment components after fusion? If two states belong 
to the same superselection sectors, they are related by 
symmetry, which protects them from decoherence. This 
may not be the case when they come from different su- 
perselection sectors, though. 

For example, consider the states |a') (g) \a^b) for i and 
j between and p — 1, where the kets will denote single 
anyons in this paragraph and the next. States of different 
j are all in the same conjugacy class, but states of dif- 
ferent i are grouped into conjugacy classes of q elements 
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(except for i = 0, which is its own conjugacy class). In 
total, we are talking about states. 

These states fuse into the states with flux for < 
k < p. The resulting states may also have one of q electric 
charges. In total, we fuse into a space containing pq 
states. Since pq < p^, what must be happening is that 
different conjugacy classes are mapped to states that are 
orthogonal in the environment subspace. 

Note that the decoherence seems to occur when fusing 
out of a state made up of different superselection sectors. 
However, fusion is the only operation that could have 
measured the relative phase between the sectors, and it 
clearly does not. Therefore, it is acceptable to assume 
that the decoherence occurs as soon as states are mapped 
into different superselection sectors. 

Returning to the question of alternative implementa- 
tions of the projection onto |0)^, it is clear that the pro- 
cedure described above does not achieve its goals without 
causing decoherence in the general case. However, in the 
special case when q = p — 1, the non-trivial powers of 
a form one conjugacy class. Therefore, the above trick 
can produce a projection onto |0)^ using only magnetic 
charges. Of course, q ^ p — 1 only holds for G ^ S3. 

For other groups, the operation could become useful if 
we could tell into which superselection sector the state 
decohered, producing a probabilistic projection onto a 
smaller space. The smaller projections may also be com- 
putationally powerful. However, since we have completed 
a universal gate-set without the results of this subsection, 
we shall work on proving universality from the previously 
constructed gates, rather than pursuing this matter fur- 
ther. 



IV. GATE-SET UNIVERSALITY 

The goal of this section is to prove the universality 
of the following qudit gate-set, which includes measure- 
ments: 

1. Controlled-X, 

2. Probabilistic projection onto |0), 

3. Probabilistic projection onto |0), 

4. Probabilistic projection onto |0)^, 

where we assume that the qudits are of dimension d > 2, 
with d prime. The first requirement on d is needed to 
make the gate-set universal, whereas the second one will 
allow us to relate this gate-set to Gottesman's gate-set 
. The above gate-set must be supplemented by a con- 
trolling computer capable of universal classical computa- 
tion. 

The above gates were selected as those arising natu- 
rally from the anyons based on the groups Zp x g Zg . The 
proof of universality of the above gate-set is the last step 



needed to show that universal quantum computation is 
feasible with these anyons. 

The proof of universality will proceed in two steps. In 
the first step we will turn the second and third gates into 
proper measurements in the Z and X bases. Most of 
the methods of the first step were described while build- 
ing computation with non-solvable anyons [l)|. The sec- 
ond step involves using the probabilistic projection onto 
|0)^ to construct magic states that complete the univer- 
sal gate-set. This is the new element needed to achieve 
universality with solvable anyons. 



A. Non-destructive measurement of Z and X 

By the end of this subsection we will have constructed 
measurements in the Z and X bases. These measure- 
ments will be non-destructive in the sense that if result 
i was obtained, the measured qudit will be in state |i) 
or \i) respectively. Because the measurements in ques- 
tion are complete, the non-destructive requirement can 
be achieved by having ancillas for every eigenstate of X 
and Z, and then using the controUed-X to swap the an- 
cillas into the computational space. 

The construction begins by producing a set of basic 
ancillas. Along the way we will also produce the X and 
Z unitary gates. 

1. \0) and \0) ancillas 

Clearly, given |0) ancillas we can use the third gate to 
produce |0) ancillas. Similarly, given |0) ancillas we can 
use the second gate to produce |0) ancillas. Therefore, if 
the initial state of the quantum computer overlaps with 
either state, we can produce both kinds of ancilla. 

Usually, the initial state of the quantum computer is 
|0). However, by using the controUed-X gate, in combi- 
nation with the projections onto |0), we can obtain these 
states no matter what the qudits are initialized to. The 
procedure is just to apply a controlled-X~^ (equivalent to 
d—1 controUed-X gates) to two qudits, and then project 
the target to the |0) space. If the initial state had some 
overlap with any of the states |j) (g) then this produces 
the desired ancillas. Furthermore, even if we allow states 
that are initially entangled, once we involve more than 
d qudits, at least one pair must have an overlap with 
the diagonal states. Therefore, |0) states can always be 
produced. 

Henceforth, we shall assume an ample supply of |0) 
and |0) ancillas. 



2. |1) states, jl) states; X gates, Z gates 

The next step is to produce |1) and |1) ancillas. The 
importance of these ancillas is that they will break the 
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symmetry currently present in the one-qudit Hilbert 
space. 

There are two symmetries in the Hilbert space that 
are not fixed by the basic four gates of our set. The first 
symmetry is a relabeling \ix) \i), calculated modulo 
d, for some < x < d. The second, is the relabeling 
Zy ^ Z, for integer < y < d. For fixed x, the second 
symmetry is a relabeling of our d*^ root of unity uj by 



Lu and a rclabcHng \jy) 



\J) 



Therefore, given an ancilla in a state \x), with x > 0, 
we can just rename it so that it becomes a |1) ancilla. 
Similarly, given an ancilla in a state \y), y > 0, we can 
relabel it as |1). In fact, both can be done simultaneously 
in a consistent fashion, even if we do not know the values 
of X and y. 

The initial states \x) and \y) can be obtained from two 
maximally mixed states. The maximally mixed states 
can be described either as a state \x) with x chosen at 
random, or a state \y) with y chosen at random. There- 
fore, two maximally mixed states serve our purpose as 
long as we do not obtain a; = or y = 0. These two bad 
cases will be detected below, in which case the process 
can be restarted with two new mixed states. 

To produce the maximally mixed states wc apply a 
controUed-X with |0) as source and |0) as target. The re- 
sult is a maximally entangled state, which can be turned 
into a maximally mixed state by discarding one of the 
two qudits. Two of these mixed states will serve as our 
ancillas. 

Given our two ancillas, which we have now labeled |1) 
and |1), we can build X and Z gates which arc consis- 
tent with the new labeling. The X gate is clearly just a 
controlled-X with a |1) state as control, whereas the Z 
gate is just a controUed-X with a |1) as target. The less 
familiar second construction is just a specific case of the 
following circuit: 









X 









\i-j}=Z^\i} 



At this point, if we were unlucky enough to get a; = 
or y = 0, then one of the transformations X or Z will 
be the identity operator. This can easily be checked by 
applying them to |0) or |0) ancillas and then using the 
available probabilistic projections. 

The X and Z gates can also be used to produce a 
reservoir of |1) and |1) ancillas that will be consistent 
with the original states. Two elements in the reservoir 
can also be compared, for example, by applying a Z built 
from one ancilla followed by a Z~^ built from the other. 
Therefore, even if the states were to decay over time, by 
using majority voting the damaged states can be weeded 
out. 

In some cases, the one-qubit Hilbert spaces do have 
natural |1) or |1) states, which implies a natural way 



of measuring or obtaining such states. For those sys- 
tems, either the natiiral ancillas or the arbitrary ones 
constructed above can be used. For example, for the 
anyons |1) = \aba~^). However, choosing a different |1) 
state is equivalent to choosing a different element a. 



3. Measurements of Z and X 

At this point all the elements are in place to produce 
measurements in either the Z basis or the X basis. 

The key element of the X basis measurement is the 
circuit 



|0) 



Z'J 



l« - i) 



applied to a |0) ancilla and the state to be measured. If 
the above circuit is repeated many times, each time with 
a different |0) ancilla, and with j varying from to rf — 1, 
we obtain the transformation 

^/3i|i) — > ^ A K) |i) • • • (g) O 1^^) • • • 

i i 

oii^ir+i). (59) 

A probabilistic projection onto |0) can then be applied 

to each qudit. If one of the qudits of the form \i — j) 
projects onto the space |0), then the outcome of the mea- 
surement is j. 

Note that because of the one-sided error model of the 
probabilistic projection, an erroneous measurement re- 
sult can never be obtained, no matter how small ppp is. 
The worst possible outcome is that after all the qudits 
have been measured, no conclusion can be reached. Of 
course, a standard small two-sided probability of error 
can also be made exponentially small by using enough 
qudits in the above measurement. 

The measurement in the Z basis proceeds similarly, 
where the transformation 

^OLi\i) — >^ai\i) ■ ■ ■ ®\i - \) ®\i - \) ■ ■ ■ 

i i 

(^\i-d+l) (60) 

is performed using the X and controlled-X gates, fol- 
lowed by a probabilistic projection onto |0). 

Finally, the above measurements can be performed 
non-destructively, by projecting all but one of the qudits. 

Alternatively, the eigenstates of X and Z can be directly 
constructed from these gates and |0) or |0) eigenstates. 
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B. Completing the gate-set 

So far, we have only shown that our gates can real- 
ize operations in the ClifFord group. In order to achieve 
universal quantum computation we need to complete the 
gate-set with an operation outside the Chfford group. 

It was shown in Ref. that the Toffoli gate, combined 
with measurements in the X and Z bases, is universal 
for quantum computation. Therefore, a successful con- 
struction of the Toffoli out of our gate-set will prove it 
universal. The Toffoli gate will be constructed out of the 
previously described operations, together with the thus 
far unused probabilistic projection onto |0)^. 

In addition to producing measurement gates, proba- 
bilistic projections are particularly useful for preparing 
magic states, which are ancillas whose use allows us to 
perform new gates such as the Toffoli. In particular, we 
shall show that we can produce the two magic states 

I0M2) = (61) 

where 5i^j is the Kronecker delta function. The first of 
these states produces the Toffoli gate up to some errors in 
the Clifford group. The second magic state allows us to 
correct these errors, and in fact, allows the construction 
of the complete ClifFord group even without the use of 
the first magic state. 

We shall begin by discussing how to use each of the 
magic states, and then afterwards turn to the task of 
describing their construction out of the available opera- 
tions. 



1. Using |<j!>Mi) 

The magic state \4>mi) and its use inproducing the 
Toffoli gate was first introduced by Shor Q , and general- 
ized to qudits in Ref. Q . We shall give a brief description 
of its use in order to give an account of the exact Clifford 
group operations needed in the last step as corrections. 

The procedure begins with a general state 

I*) = ^ V'a,fc,c|a)<8) |fo)<8) |c), (62) 

to which an ancilla \(i>Mi) is appended. A controUed-X^^ 
is applied to the first data qudit with the first ancilla qu- 
dit as control. Similarly, a controlled-X~^ is applied to 
the second data qudit from the second ancilla qudit, and 
a controUed-X is applied to the third ancilla qudit, from 
the third data qudit. The first two data qudits are then 
measured in the Z basis and the third data qudit is mea- 
sured in the X basis. If the results of the measurements 



are a, (3 and 7 respectively, then the remaining qudits 
are left in the state 

i'aAc^'" \a-a)®\b- I3)®\{a- a){b -(3) + c). 

a,b,c 

(63) 

The corrections begin by applying a.n X°'(^X^(^X 
gate followed by a controUed-X^ from the first qudit to 
the third qudit and a controUed-X" from the second qu- 
dit to the third qudit. The state then becomes 

Y,i'a,b,cUj'"'\a)(g)\b)(E)\ab + c). (64) 

a,b,c 

All that is needed to complete the Toffoli gate is a Z~'^ 
gate applied to the third qudit and a phase uj''"-'' applied 
to the first two qudits. Unfortunately, we must first build 
the latter transformation out of the second magic state. 

2. Using \<pM2) 

Once again, the magic state is appended to a pair of 
qudits. Now controUed-X gates are applied with the data 
qudits as source and the ancilla qudits as targets. Then 
the ancilla qudits are measured in the computational ba- 
sis. The outcomes a and /3 will be uniformly distributed, 
and at the end we will have produced the transformation 

^V'a,6|a)® |6) -^^V'a.bW^"-*"-^ \a)®\b). (65) 

a, 6 a,b 

This procedure randomly and uniformly chooses a 
computational basis state and multiplies it by a phase 
of uj. Repeated application of this transformation will 
eventually yield any of the d"^ states of the form 

^Va,bC^^('''')|a)0|6), (66) 
0,6 

where / is an arbitrary integer-valued function. This 
process is effectively a classical random walk on a di- 
mensional periodic lattice with d nodes, where each use 
of a magic state is equivalent to taking one step. Because 
the lattice is finite, after a polynomially large number of 
steps the probability of not having arrived at least once at 
any one of the above states becomes exponentially small. 

The final correction needed to complete the Toffoli 
gate was the phase transformation to the state with 
/(a, b) — jab, and can therefore be realized using many 
copies of the second magic state. All that remains to 
prove universality is to describe the production of the 
magic states. 

3. Making the magic states 

The final piece of the puzzle is the production of the 
magic states using the probabilistic projection onto |0)^. 
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Probabilistic projections onto a subspace are particu- 
larly powerful for making magic states, because it can be 
assumed that they successfully project into the subspace 
every time. That is, if the probabilistic projection does 
not project onto the desired subspace, the state is tossed 
out and the procedure is restarted from the beginning. 
Therefore, the probabilistic projection onto lO)""" effec- 
tively takes a state and removes the |0) component of 
the state: 



d-1 



d-1 



^ai\i) — > A^ai\i}, 



(67) 



1=0 



where A is some normalization constant. In fact, by com- 
bining this projection with the X gate, we can remove 
any of the components 

The main strategy for this section is to construct a se- 
ries of ancilla states of increasing complexity, until finally 
the desired magic states are obtained. At this point, we 
have a supply of ancillas of the form \i) and \ j) for any 
i and j. From the |0) state we can also make the an- 
cilla (|0) -I- |l))/\/2 by removing all \i) for i > 1 with the 
probabilistic projection. 

The next step is to produce entangled two-qudit ancil- 
las. Given a supply of ancillas of the form \^) = J2i V'il^) 
we shall produce ancillas of the form 

/ d-1 \ 



i=l 



I*') = (^^o|0)®|l) + ^^i|^)® |0) 

d-1 



(68) 



i=0 



The procedure begins with the state 



I*) ® (^|0) + ^|1)) = ^ EEV'iK) ® (69) 

which in general has 2d non-zero coefficients. We need 

to remove d of these coefficients to obtain the state 

The procedure, done once for each k from 1 to d — 
1, is the following: First, apply a controUed-X*^ with 
the left qudit as source and the right qudit as target. 
Then, the right qudit is projected onto |0)^, and finally 
the controlled-X*^ is undone. For each fc, we remove the 
components |0) ® |0) and \ — l/k) ^ |1). The operation 
— 1/k is modulo d as usual, and ranges over all integers 
between 1 and d—1 because d, is prime. Therefore, given 
a supply of l^*) ancillas, we can probabilistically convert 
some of them into a supply of \^') ancillas. 

Note that the above procedure works even if the coef- 
ficients ipi represent the state of other qudits, as long as 
these are ancilla qudits that can be tossed out if the pro- 
jection procedure fails. In the same spirit, given ancillas 
of the form 



(70) 



we can produce the three-qudit ancillas 

1 1 

I*') = EE<^^'^-|^) ® b') ® (71) 
i=o i=o 

The procedure again involves appending (|0) -f |l))/-\/2 
to the ancilla |$), which now generically has eight non- 
zero coefficients, and removing four of them. This is done 
with a set of controUed-X gates with the third qudit 
as target, followed by a probabilistic projection of the 
third qudit onto |0)^, followed by the inverse controUed- 
X gates. If we use two controlled-X^^ gates controlled 
by the first two qudits respectively, the projection will 
remove the components with labels |0)|0)|0), |l)|0)ll) 
and |0)|1)|1). In addition, using two controlled-X^'^"^-'/^ 
gates, we remove |1)|1)|1) and |0)|0)|0) (again). These 
are the four states that need to be removed to produce 
the ancilla |$'). 

The above two procedures allow us to finally produce 
the desired magic states. Starting with |0)(g) |0), we apply 
the first procedure to each ancilla and then apply the 
second procedure to the appended qudits. The resulting 
state is 

^ d-1 d-1 

E E I*) ® li) ® \Si,o) \6j,o) ® \Si,oSj,o). (72) 

i=0 j=0 

If the last three qudits are measured in the X basis, and 
the results are 0, and 1 respectively, then we will have 
produced the magic state \(j)M2)- 

In fact, measuring in the X basis and only accepting 
if the result is zero is a convenient way to unentangle 
the system with temporary qudits. Therefore, the pre- 
viously described procedures can be combined into the 
probabilistic transformation 



d-1 d-1 

EE^^'j- 

i=0 j=0 



i) O |j) 



d-1 d-1 

EEV^^: 

i=0 j=0 



(73) 

where the first state is either transformed into the second 

state with some nonzero probability, or else it is damaged. 
The above transformation has only been discussed so far 
for n = m = 0, but a trivial use of X gates before and 
after the transformation will allow any n and to. 

Starting with |0) (8) |0) (8) |0), repeated application of the 
above procedure can produce 

jEEN)®li)®|0) 



i=0 j=0 



-1 d-1 



3EEi^)®ij')®|o)(8)(8)i'^^."'^j-.™) 



d 



i=0 j=0 n=0 m=0 

d-ld-1 d-1 d-1 

;y E E 1^) ® \3) ® \^i,nhm), (74) 



=0 j=0 



d 



i=0 j=0 



n=Q m=0 
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where the second step involves only controlled-X gates 
from the extra qudits to the third qudit. Erasing the 
extra qudits with a measurement in the X basis, and 
retaining only when all results are zero, produces the 
desired magic state |0mi)- 

The construction of the magic states out of the prob- 
abilistic projection onto |0)^ completes the description 
of the Toffoli gate. Though the procedures of this sec- 
tion are far from optimal in terms of resources, they are 
sufficient to demonstrate universality. In particular, this 
completes the proof that universal quantum computation 
is feasible with anyons from groups of the form Zp XgZg. 



V. COMPUTATIONAL POWER OF MAGNETIC 
CHARGES 

In this section, we will be interested in classifying the 
computational power that can be achieved by braiding 
anyonic magnetic charges of a finite group. The range 
of operations that can be achieved by braiding is closely 
related to the structure of the group to which the mag- 
netic charges belong. In particular, the possibility of re- 
alizing the operations of controUed-X and Toffoli (cquiv- 
alently a doubly-controUed-X) are respectively related to 
the group properties of nilpotency and solvability. These 
standard properties of group theory will also be defined 
below. 

There are certain important assumptions that go into 
the discussion in this section. First, we assume that each 
qubit is carried by a pair of anyons. Furthermore, we 
choose a computational basis corresponding to the states 
of definite flux (e.g., |0) = \g) for some g G G). We 
remind the reader at this point that the state l^r) corre- 
sponds to an anyon of magnetic charge g paired with a 
compensating anyon of charge g^^ whose only purpose 
is to allow the pair to move through the system without 
introducing undesired correlations. Finally, we will re- 
strict the discussion to operations that can be achieved 
by braiding magnetic charges. The consequences of lift- 
ing these restrictions will be discussed near the end of 
this section. 

Let the fluxes corresponding to the zero and one states 
be the elements b, b' G G respectively. If we desire a 
coherent superposition between the zero and one states, 
they must be in the same conjugacy class, and therefore 
b' = aba~^ for some non- trivial a € G. This is summa- 
rized by: 



The most general conjugation possible is by a function of 

the form 

fid) = Cigc2gc3g ■ ■ ■ gcn (76) 
= (rfi.g'dri) {d^g'd^') ■ ■ ■ (d„_ig'd;ii) dn, 

for some fixed elements {cj} € G. In the second line, the 
expression has been rewritten in terms of g' = gb^^, and 
new elements {di} G G which can easily be determined 
in terms of {cj} G G. For example, ^2 = Cibc2. 

The power of the second line is that it expresses the 
conjugation as a composition of two basic operations. 
The first is a conjugation by an ancilla with flux dn, and 
is independent of the state of the first qubit. The second 
is conjugation by a product of conjugates of g' , which was 
defined so that if g = b then g' = 1 and the product of 
its conjugates is trivial. In the other case, if g = b' then 
g' = [a,b] = aba~^b~^, and the operation is conjugation 
by a product of conjugates of [a, b]. 

We define Cg{x) as the conjugacy class of x in G, and 
Cq{x) as the group generated by the elements in Cg{x). 
The operations discussed so far are conjugation by fixed 
elements in G, and controlled conjugation by elements in 

e*{[a,b\). 

The most natural controlled operation is the logical 
controUed-X gate, which acts as a controlled conjugation 
by a. Naturally, if ^ 1, then we could arrive at the 
qudit state |2) = |a^6a~^). However, our interest lies in 
proving that certain groups cannot produce a controUed- 
X, in which case it is sufficient to prove that a controlled 
conjugation by a is unfeasible. 

It seems that a requirement for a controlled conju- 
gation by a is the existence of elements a, 6 such that 
a G ([a, b]). There is a potential loophole in the argu- 
ment, though, because different qubits could use different 
basis fluxes. The target qubit could use 62 as the zero 
state, and a2&2a2^^ as the one state. If 02 G CQ([ai,6i]) 
then the controUed-X would be possible. Considering 
many qubits requires a sequence of non-trivial elements 
{ai} and {bi} which satisfy, at a minimum, the condi- 
tions: 



ttj+i G C^{[ai,bi])- 



(77) 



The above equations are related to the series of sub- 
groups of G, deflned by 



(5(0+1)) ^ [g!(0)),G 



(78) 



|0) = 



ID 



\aba 



(75) 



Even if the basis in use is a qudit basis, with additional 
states, we will only concern ourselves with states that 
have support on the above two basis vectors. 

Consider now a pair of these states. We are interested 
in the operations that can be achieved by conjugating the 
second state by the flux of the first state, with the help 
of ancillas. Let g G {b, b'} be the flux of the flrst state. 



with base case G((°» = G. By deflnition, if G G^^-'" 
then [a„6,] G G^^^+i^. Furthermore, since the group 
q{{3+'^-)) is normal in G, the requirement on Oj+i reads 



(79) 



Of course, oi G G^^°^^ = G. Therefore, repeating the 
above argument shows that a controlled-X requires a, G 
G((*-i)) with a^^l for every i > 1. 
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Given that G is finite, and G((^+i» C G«^», the series 
must converge after a finite number of subgroups to some 
final subgroup 

G((oo))_ Tj^g fjjjal sub group can either be 
trivial or non-trivial. The groups with G((°°» = {1} are 
called nilpotcnt. The conclusion thus far is that nilpo- 
tent groups cannot implement a controUed-X by braid- 
ing. The inverse of this statement — i.e., that groups that 
are not nilpotent can implement the controUed-X — will 
be shown in Section IVll 



A. Conjugations with multiple sources 

A similar analysis can be used to study the relationship 
between group structure and gates produced using mul- 
tiple qubits as sources of conjugation. Clearly any group 
that is not nilpotent can produce a series of controlled- A" 
gates with different sources. However, certain groups are 
capable of producing much more powerful gates such as 
the Toffoli, which is universal for classical computation. 

In the rest of this section we shall prove that groups 
that are solvable cannot produce a Toffoli gate, or equiv- 
alently universal classical computation, by braiding mag- 
netic charges. This connection between universality for 
classical computation and non-solvability had been pre- 
viously identified by Barrington ^3 i'^ 1989. Though we 
shall mostly be interested in groups that are solvable, this 
result will place limits on the power that we can expect 
to obtain from braiding magnetic charges. 

Just as above, the most general conjugation with m 
sources is the conjugation by a function of the form 



■ ■ ■ (d^^,gl_^d--i) dn, (80) 

where g[ — gib^^ and the indices i take values from 1 to 
m. For brevity, we assume that all qubits are expressed 
in the same basis, though the general case would not be 
very different. 

The Toffoli gate is simply a conjugation by a function 
/T(ff 1,52)1 such that 



/T(c^c') 



kl 



which has been expressed as a function of and where we 
introduced c = [a,b]. In order to produce the Toffoli gate 
using conjugation alone, we must be able to express the 
above equation in the form of Eq. (|80|l with m = 2. We 
shall show that this is not possible for a solvable group. 

For m = 2, Eq. (|80|l is a product of conjugates of 
g'l and 52- We can rewrite it by moving all the conju- 
gates of g'l to the left, and all the conjugates of (73 to 
the right. In the center we will pick up factors of the 
form [dig'idj^ Tdjg'2d~^] and commutators of commuta- 
tors, and so on. In the end we will obtain: 



where dn is a constant element of G, fi is a product of 
conjugates oi g[, and /c(5i,52) the factor with all the 
commutators. The function fc has the property that 

/c(5i,l) = /c(l,ff^) = l. 

Setting f ~ fr implies the conditions 

1 = /(l,l)-d„, 

1 = /(c,l) = /i(c)d„, 

1 = /(l,c)-/2(c)d„, 

a - /(c,c) = /i(c)/c(c,c)/2(c)d„, (83) 

which imply /c(c, c) — a. 

However, fc has the additional property that, if N is 
a normal subgroup of G containing c, then fc{c,c) G 
[N, N]. Furthermore, since c = [a, 6], the requirement on 
c needed to express the Toffoli function in product form 
is 



ce[N,N], 



(84) 



for any normal subgroup N . This condition is related to 
the series of subgroups defined by 



(85) 



/(3l,52)=/i(5i)/c(5i,52)/2(52K 



(82) 



again with base case G'-^-' = G. Just as before, this series 
must converge to a final subgroup G(°°). The groups 
where G^°°^ = {1} are known as solvable. Any group 
that is nilpotent is also solvable. 

Because the subgroups are all normal in G, the 
requirement of Eq. H84I) can only be satisfied if c, which 
by definition cannot be 1, is contained in G'°°'. We have 
therefore shown that if the group is solvable, then the 
function fx cannot be expressed in product form, and 
therefore we cannot conjugate by it. This is true even 
if the target of conjugation is in a known state, which 
implies that even if we had used the target as a source 
of conjugations as well (i.e., by using it to conjugate an- 
cillas, and then using the ancillas) the Toffoli gate would 
still not be feasible by using only braiding of anyons from 
a solvable group. 

The fact that the Toffoli gate can be produced for 
non-solvable groups is a consequence of the results of 
Refs. and will not be discussed here. In fact, the 

computational model discussed in this section resembles 
the non-uniform deterministic finite automata presented 
in Ref. TT*!. For non-solvable groups, the two models are 
almost identical. Nonetheless, for solvable groups, the 
magnetic charges presented in this section have signifi- 
cantly less computational power, because the zero and 
one states have to be represented by group elements in 
the same conjugacy class. 



B. Summary of computational power 

The results discussed so far have been summarized by 
Tabled For each type of group, it describes the computa- 
tional operations that can be achieved through braiding 
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Abelian 


Nilpotent 


Solvable 


Example 


Computational Power 


yes 


yes 


yes 




I 


no 


yes 


yes 


Q 


X 


no 


no 


yes 


Sa 


ControUed-X 


no 


no 


no 


As 


Toffoli 



TABLE I: Computational power achieved by conjugation for 
different groups. 



of magnetic fluxes, as well as an example. The examples 
are the smallest group in the class, with the exception 
of the abelian case where the trivial group could also be 
listed. For the non-abelian, nilpotent case there are two 
examples with eight elements: the dihedral group D4, 
and the quaternionic group Q, which is listed in the ta- 
ble and has elements ±1, ±i, ±j, ±fc. 

The most basic case is when G is abelian, in which case 
it is also nilpotent and solvable. Clearly conjugation can 
only produce the identity transformation. In fact, ev- 
ery superselection sector consists of a one dimensional 
Hilbert space, and therefore quantum information can- 
not even be stored in abelian anyons in a topologically 
protected manner. 

At the other extreme are anyons from non-solvable 
groups. Universal classical computation can be accom- 
plished through braiding, and universal quantum compu- 
tation can be obtained by completing the gate-set with 
measurements in the X and Z bases. The complete con- 
struction for this case is described in Ref. ij. 

Anyons from groups that are solvable but not nilpo- 
tent, can also be used for universal quantum compu- 
tation, but the construction is more complicated. A 
controUed-A' can be constructed from flux braiding, and 
measurements in the X and Z bases can be constructed 
in a manner similar to the non-solvable groups. How- 
ever, to complete a universal gate-set, fusions of electric 
charges must be employed. The proof of universality, 
along with the details of the gates, will be the subject of 
the rest of this paper. 

Finally, anyons from groups that are nilpotent seem 
insufficient for universal computation. In the construc- 
tions for the non-nilpotent groups, the only operation 
that can produce entanglement between multiple qudits 
is the controUed-A or Toffoli gates obtained by braid- 
ing fluxes. However, for nilpotent groups, braiding fluxes 
does not seem to yield an operation capable of produc- 
ing entanglement. Either a new type of operation, or a 
different basis must be used. Simple modifications to the 
basis, such as encoding a qudit on multiple anyons, are 
of no help. However, there are countless strange bases 
that are hard to discredit. For example, a lattice of elec- 
tric charges could serve as a Hilbert space, with magnetic 
charges used to create or measure entanglement among 
the charges. Therefore, while the prospects of univer- 
sal computation with nilpotent anyons seem bleak, the 
question remains open. 



VI. SOLVABLE NON-NILPOTENT GROUPS 

In this section, we will prove that anyons based on a fi- 
nite group that is solvable but not nilpotent are sufficient 
for universal quantum computation. The first step will 
be to decompose an arbitrary group G that is solvable 
but non-nilpotent, into a form similar to the previously 
studied Zp XgZq groups. The proof of universality will 
then be a small generalization of the ideas presented in 
Section Iml 



A. Group decomposition 

Let G be a group as above and define H = G^'^°°^^ in 
terms of the series discussed in Section [V] Because G is 
non-nilpotent H is non-trivial, and because G is solvable 
H ^ G. Furthermore, H is normal in G and G/H is 
nilpotent. The second fact is due to 



(G/H) 



((«+!)) 



iG/Hf''>\{G/H) 

G^W^gI /i? = G«*+i»/i?, (86) 



and therefore (G/iJ)((°°» = H/H = {1}. 

Any nilpotent group can be written as the direct prod- 
uct of its Sylow p-groups, which are groups whose order 
is a prime power. Therefore 



G/H = Kg, X Kg. 



X K, 



(87) 



where Kq denotes a group of order q^' 
and integer to. We further define NKq 
Kq^ to the fuU group G that is NKq./N 



for some prime q 
to be the lifting of 
Kq^. Note that 



to maintain consistency with the notation in Section [llll 
the primes involved in these p-groups are labeled by the 
letter q. 

Having fully characterized G/H , we turn to the study 
of H itself. Let N be the largest normal subgroup of G 
that also satisfies N C H and N ^ H . If more than one 
subgroup satisfies the above requirements, then let N be 
any such subgroup. Because H is finite, there must be 
at least one maximal subgroup. 

We shall prove that H/N — for some prime p and 
integer n. The basic idea is that working modulo , H/N 
is a normal subgroup of G/N . Furthermore H/N has no 
proper subgroups that are normal in G/N . In particular, 
this implies that H/N is abelian, because its commuta- 
tor subgroup is a normal subgroup of G/N. Note that 
the possibility that the commutator subgroup of H/N is 
equal to H/N is excluded because H/N is solvable. 

For any x G H/N consider C^^jy(x), the group gener- 
ated by the conjugates of x in G/N. This is a subgroup 
oi H/N and is normal in G/N. Therefore 



C*^j^ix) = H/N, Vx e H/N, 



(88) 



which implies that all elements in H/N, with the excep- 
tion of the identity, have the same order. That is because 
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conjugates of x have the same order as and a product 
of elements of order fc, in an abehan group, must have 
order less than or equal to k. This concludes the proof 
that H/N = 17^. 

Thus far, we have the following tower of groups 



N CH C HKq^ C G, 



(89) 



where N ^ H and HKq. are all normal in G, and the group 
HKq- can be any of the groups found above. 

Because (G/7V)»°°)) = H/N = Z^J, the group G/N is 
also solvable and non-nilpotent. However, its structure 
is simpler than that of the full group G. We shall there- 
fore be interested in working modulo N , and shall denote 
groups modulo by a tilde. That is 

G = G/N, HKq^ = HKqjN, H = H/N = Z]]. (90) 

The final step is to study the relationship between 
H and the groups HKq. . By construction, we know 
[G,H] = H = Z^,hut what about [HKq^,H]7 Because 
both HKq. (for any i) and H are normal in G, [HKg^ , H] 
is normal in G and, furthermore, is contained in H. But 
N was defined to be the largest proper subgroup of H 
that was normal in G. Therefore H has no proper sub- 
groups that are normal in G, and [HKq. , H] must be 
either the trivial group or all of H. 

li Qi = p then HKq^ is a p-group, and therefore nilpo- 
tent. This means that [HKp, H] ^ H and by the previous 
paragraph [HKp, H] = {1}. The rest of the groups HKq^ 
can either commute or not with H. However, because 
[G, H] — H, at least one of them must not commute. 
Fix an i such that [HKq. ,H] — H, and define K = Kq. , 
HK = HKq. and q ~ qi. This will be the group to take 
the place of Zq. 

We would like to show that there exists an element 
b e HK, such that [b, H] = H. Let X be the stabilizer 
of H in HK, that is, the largest subgroup of HK such 
that [X, H]^l. Clearly H C X and X HK. Because 
HK/X is nilpotent, it has a non-trivial center. Let b G 
HK be any element that projects, modulo X to one of 
the non-trivial elements in the center. We will show that 
[b, H] is normal in G, which implies [b, H] = H. The 
proof is that modulo X (which is normal in G), every 
element g (z G commutes with b. Therefore gbg^^ — bx 
for some x G X and 



9 [b. h] g' 



bxh'x-^b-^h'~^ 



bh'b-^h'-^ e b,H 



(91) 



for any h G H, where h' — ghg^^ G H. 

To summarize, working modulo N, we have the follow- 
ing tower of subgroups: 



H C HK C G, 



(92) 



with H = Z'p for some prime p. Furthermore, HK/ H — 
isT is a subgroup of order a power of q, for some prime q 
not equal to p. Finally, 3b G HK such that [b, H] = H. 



Note that this notation is consistent with the one used 
in SectionUm That is, if G = Zp Xg Z,, then TV = {1}, 
H = Zp, K = Zq, and the definitions of p, q, and b would 
be consistent. 



B. Examples 

There are a few good examples to keep in mind that 
illustrate the potential new complications arising from 
groups with more structure than Zp xg Zq. 

The first example is A4 = Z| Xe Z3. The group can 
be described as oi = (12)(34), 02 = (23)(41), b = (123) 
with 



1, 
1, 



aia2 — 0201, 



ba\a2b 



i+3 



a, a 



l"'2 



(93) 



For this group N ^ {I}, H = Zl and K = Z3. Its 
most important feature is that p ~ 2, which was not 
previously possible. Because p — 2 implies working with 
qubits, these groups will be have to be handled specially. 

The next example is G = (Z3) Xg (Z3 x Z2). Let ai, 02 
be the generators of Z§, let b be the generator of Z2 and 
x be the generator of the remaining Z3. The semidirect 
product is defined by the conjugations 



''1^2'- 



^2 1 



(94) 



For this group H 



Z§ because [G,G] 



[G,H] 



H. 

The subgroup generated by aia^^ is normal in G and 
therefore = Z3. Finally H/N = Z3 and K = Z2. Note 
that X commutes with H modulo N, as discussed in the 
last section. 

The final pair of examples illustrate the case where K 
is non-abelian. The examples are Z§ xg Q and Z| Xg D4. 
Labeling the generators of Z| by ai and 02, the semidirect 
product for Z3 x ^ Q is defined by 



n"2 



1 _ „^+i/„^-s/ 



a 



On 



(95) 



where ±1, ±i, ±j, ±fc are the standard quaternionic el- 
ements. For Z| Xe -D4 the semidirect product is defined 
by 



f3alalP 1 = a\a2 ^, 7<a27 ^ = a\a^. 



(96) 



where the relations = 7^ = 1 and 7/^7 — /3 ^ define 
D4. 

In both of the above cases p = 3, q = 2, N = {1} and 
H — Z\. However, for Z\xg Q the non-trivial elements 
of H are conjugate to one another, and none of the non- 
trivial elements of Q commute with any of the non-trivial 
elements of H. The Z\ Xg D4 case divides H into three 
conjugacy classes (including the identity). Furthermore, 
each of the elements of the form /3*7 commute with two 
non-trivial elements of H. These differences will become 
important when discussing the operations involving elec- 
tric charges. 
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C. A'^-invariant ancillas 

The first lesson from the above analysis is that we 
should work modulo N. That is, we want flux states la- 
beled by elements oi G = G/N, that are invariant under 
N. The idea of iV-invariant states was already discussed 
m Ref. when generalizing simple non-abelian anyons 
to non-solvable ones, and therefore the discussion below 
will be brief. 

A basis for the iV-invariant magnetic fluxes is just \g) 
for g ^ G. The braiding and fusion properties of these 
states behave almost exactly as if the full group were 
G and these states were flux eigenstates. The only dif- 
ference is that when fusing two anyons from pairs with 
opposite fluxes, the probability of disappearing into the 
vacuum is lower. 

Even producing anyons from the vacuum behaves cor- 
rectly with respect to A^-invariance. Pairs produced from 
the vacuum are naturally invariant under the full group 
G. Normally, when braiding with other states, this in- 
variance will be broken. However, if the vacuum pair 
only interacts with A'^-invariant states, then the invari- 
ance under the group A^ will remain. 

At this point we will change our requirements for the 
physical system. Instead of requiring a reservoir of flux 
ancillas for every element of G, we will require a reservoir 
of A^-invariant flux ancillas for every element of G. This 
is likely a reasonable modification, as it appears that the 
latter ancillas are no harder to produce than the original 
ones. 

It should be noted that, when working modulo A^, the 
electric charges need no modification. That is, because A^ 
is normal in G, any representation of G extends to a rep- 
resentation of G that is invariant under A^. Furthermore, 
fusing two A"-invariant electric charges must produce a 
new A'^-invariant electric charge. Therefore, working with 
A'-invariant electric charges simply involves working with 
a subset of the charges of the group G. 

Given the above caveats, we can effectively replace the 
group G with the group G = G/N, which will be done 
without further comment for the rest of this section. 



D. Computational basis 

We will begin by defining an extended computational 
basis, and discuss the operations that can be performed 
on this extended subspace. Toward the end of this sec- 
tion, a subset of these states will be singled out as the 
true computational basis. 

Let oi, . . . , a„ be a set of generators for H — Z^, and 

recall the definition of the element b E G. The extended 
computational basis consists of the states 



(97) 



To prove that the states are all distinct consider the 
map from H H defined by 



i5, 



(98) 



Because H is abelian, this map is an homomorphism for 
any g G. In particular, since [6, H] — H , the homomor- 
phism defined by [b, ■] is surjective and has trivial kernel. 
That is, no element of H commutes with b. But 



hbh-^ = h'bh' 



{h'-^h)b{h'-^h)-^ = b, (99) 



for any elements /i, h' G H , which can only be true if 
h = h'. 



E. Basic operations 



The generalized controUed-X is the transformation 



1*1, ■ 



• • , in) 



(100) 



It can be implemented as a conjugation of the second 
anyon by a function of the flux of the first anyon such 
that 



f{hbh-^) = h, 



(101) 



for any h E H. Because the map [b, ■] defined above is 
just a permutation of the elements of H, it has a finite 
period (say I). The desired function is 



[gb-\b],b 



(102) 



which consists of / — 1 nested commutators. The final 
commutator needed to complete the period is the one 
formed in the expression gb~^ when g has the form hbh~^. 

At this point, one may wonder how does working mod- 
ulo a normal subgroup A^ affect the discussion regard- 
ing the computability of the controUed-X gate. The 
controUed-X can only be implemented because G is non- 
nilpotcnt. In a sense, G was constructed to be as small 
as possible, but still maintain the property of being non- 
nilpotent. On the other hand, if a group G is nilpotcnt 
to begin with, then any subgroup or quotient group will 
also be nilpotcnt, and no controlled- A can be constructed 
using braiding. 

Using the same techniques as in Section IIII Dl anyon 
fusions can be used to perform measurements. Fusion 
with \b^^) ancillas produces a probabilistic projection 
onto |0, . . . , 0). Fusing the two anyons that form a qudit 
is a probabilistic projection onto 



p-i p-i 

i,=0 i„=0 



.,«„). (103) 



where each of the i's takes values from to p — 1. 



As usual, to complete the probabilistic projection, 
these fusions must be supplemented by a reservoir of 



21 



|0, . . . , 0) and jO, . . . , 0) ancillas. The first case is trivial, 
because the existence of these ancillas has been assumed 
as one of the physical requirements of the system. The 
production of |0, . . . , 0) ancillas is more complicated and 
will occupy the rest of the subsection. 

The procedure to distill |0, . . . , 0) states begins with a 
pair created from the vacuum and a |0, . . . , 0) ancilla: 



|Vac)(g)|0, ...,0). 



(104) 



Using only braiding, an incomplete swap is applied to the 
state: 



|Vac) 



|0,...,0) 



X- 



X 



Once again, the circuit denotes the action of the conju- 
gations on the computational basis, but their extension 
to the full Hilbert space needs to be discussed. After ap- 
plying the necessary braidings to perform the circuit, the 
top state is fused with a \h~^) ancilla. If the fusion does 
not produce the vacuum state, the final product is dis- 
carded and the procedure restarted from the beginning. 
Since conjugations cannot change the superselection sec- 
tor, the only case that needs to be considered is when the 
vacuum state is created in the superselection sector that 
contains the computational subspace (i.e., the conjugacy 
class of h) . In this superselection sector the vacuum state 
has the form 



|Vac) cx |0, ...,6) + |*^>, 



(105) 



where l^""^) is a state in the space spanned by vectors of 
the form \ghg~^) that are not contained in the computa- 
tional basis. 

Because we want to guarantee that after the controUed- 
X the state |0, . . . , 0) remains in the computational sub- 
space, we need the conjugation function to satisfy 



/(G) e ff, and 
f{hbh-^) ^ h MheH. 



(106) 



The second requirement can be satisfied by choosing / 
as a sequence of commutators as in Eq. (|102f) . as long 
as the number of commutators is one minus a multiple 
of I (the period of [6, •]). Furthermore, the result after i 
commutators must be contained in G^'^*^^. Because the 
series is finite, G*^'--'^^ = G^*-""-*^ = H for some finite j, 
and the first requirement can also be satisfied by defining 
/ to be a long enough sequence of commutators. Both 
requirements can be satisfied simultaneously by correctly 
choosing the number of commutators in the expression, 
and this completes the definition of the first controUed- 
X. 



The second controlled-X can be a regular controUed- 
X because in this case the control is known to be in the 
computational subspace. In the end, the vacuum state 
will be conjugated by an element of iJ, and therefore can 
only have fiux 6 if it was originally in the computational 
subspace. 

Having completed the construction of the |0, . . . , 0) an- 
cillas, all that is required to complete a universal set of 
gates is an analog of the probabilistic projection onto 
|0)^ constructed out of fusions of electric charges. 



F. Using electric charges 

The ideal goal for this section would be the construc- 
tion of the probabilistic projection onto |0, . . . , 0)^ gate. 
Unfortunately, this is not possible for most groups. How- 
ever, we will produce a pair of gates that have an equiv- 
alent computational power. 

The first gate involves a non-trivial subgroup A d H, 
to be defined later, which could equal all of H. Note that 
this subgroup defines a subspace of the computational 
space spanned by 



\XbX- 



(107) 



for all elements A G A, which will also be denoted by A. 
The probabilistic projection onto A will be the first gate. 

The second gate is the probabilistic projection onto 
|0, . . . , 0)^ n A. This second gate can be though of as an 
application of the first gate, followed by a probabilistic 
projection onto lO)""" that only works on states contained 
in A. For the moment, we will assume that the first gate 
can be implemented, and work on the construction of the 
second gate. 

The basic building block for this section involves work- 
ing with the state to be measured \^) and an electric 
charge pair in the vacuum state |i?(/))^ of some non- 
abelian representation R. The state to be measured is 
contained in the computational basis and can therefore 
be expanded as 



I*) = ^ ^h\hbh- 



(108) 



heH 



where, as in Section IIII Gl the coefficients {iph} could 
be numbers or could denote the state of the rest of the 
system. 

Using braiding, the state j^*) can be entangled with 
the electric charges. In particular, if (l){g) is a function 
constructed as a product of g and fixed elements of G, 
then the following transformation can be realized: 



® \R{I))j^ J2 ^h\hbh'^) ® \R{^{h)))j^. (109) 



heH 



Note that the state of the electric charge can depend 
on (f>{h) rather than (j)(hbh~^) by composing with the 
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function defined in Eq. ifTU^ . That is (l){f {hbh''^)) = 

Now the electric charge pair is fused together, and the 
resulting particle is measured. More specifically, in ac- 
cordance with the discussion in Section 1111 HI we just 
check whether the resulting particle belongs to some one- 
dimensional representation labeled 7. If the charge 7 
is detected, then the electric charge will have disentan- 
gled with the state being measured, because its internal 
Hilbert space is one dimensional. Furthermore, because 
each one-dimensional representation occurs only once in 
the decomposition of R ® R* , the state will be unen- 
tangled with the environment as well. The proof of the 
latter property uses Schur's lemma and the fact that if 
|Mi)jj and \M2) always fuse into representation 7 then 
\MiM\)j^ will always fuse into the vacuum. 

The result of the complete operation, when the out- 
come 7 is obtained, is the transformation 

l*> F^{h)^-ii^h\hbh-^), (110) 

heii 

where the state after the measurement has been left un- 
normalized. The coefficients Fh^^ depend implicitly on 
the original representation i?, and will be defined care- 
fully below. 

The above procedure can be repeated many times for 
different functions (j){g). If on each occurrence the out- 
come 7 is obtained, the resulting (unnormalized) state 
will be 

J2\YlF^^,,^^Ai,f,\hbh-'), (111) 

heH / 

where $ is the set of functions used. As usual, if the 
outcome 7 is not obtained on each instance, the state 
is discarded, and the probabilistic projection reports a 
projection onto the complement. 

We assume that all functions in the set $ are products 
of conjugates of the input, and therefore (j>{I) = I for 
any S $. Because \R{T))j^ is the vacuum state, it will 
always fuse back into the vacuum. Therefore, if 7 is a 
non-trivial representation then Fj^-y = and the above 
operation projects out the |0, . . . , 0) state. 

At this point we have almost constructed a probabilis- 
tic projection onto |0, . . . , O)"*" n A. The states outside 
of A can be removed using the probabilistic projection 
onto A, which for the moment we assume can be imple- 
mented. Therefore, the desired gate will be complete if 
the coefficients 

n ^^(A)-7 (112) 

ct>e<s> 

are non-zero and equal for every non-trivial A G A. The 
requirement of equality is accomplished if the orbits un- 
der the functions in $, of all non-trivial A e A, are equal. 



More specifically, let $ be a set of maps from A to A 
that fix the identity. We say that $ is balanced on A if 
it satisfies the relation 

#(Ai ^ A') = #(A2 ^ A') V Ai, A2, A' e A - {/}, 

(113) 

where #(A A') denotes the number of elements 
such that (j>{X) = A'. The requirement that $ be balanced 
guarantees that the expressions in Eq. (|112|l are equal for 
every A. Of course, for the coefficients to be non-zero, we 
must prove separately that the value of Fx^-y is non-zero 
for every non-trivial A G A. 

The goal for the rest of this section is, therefore, to 
find: a subgroup A of H, an irreducible representation R 
of G, and a one-dimensional representation 7 of G such 
that: 

1. The probabilistic projection onto A can be imple- 
mented. 

2. F\^^ ^ for every non-trivial A G A. 

3. There exists a set of maps from A to A, that is 
balanced on A, and can be expressed as 

0(5) = n5»ff5r\ (114) 
for some elements {17^} g G. 

1. Choosing A 

There are groups, such as Z| Xq D4, for which there 
is no choice of R and non-trivial 7 such that F^^^ 7^ 
for all non-trivial h d H. It is therefore advantageous 
to choose A as small as possible. Furthermore, a small 
A will also help when proving the existence of a set of 
functions balanced on A. 

Let a be a non-trivial element of H, and consider the 
set of functions of the form 

H9)^l[9^99^\ (115) 

i 

such that (j>{a) — I. The kernel of each of these functions 
is a subgroup of H that contains the element a. We define 
A as the intersection of all these kernels. 

Because there are a finite set of maps from H to H, 
we can find a finite set of functions {(pi}, in the form of 
Eq. I|115|) . satisfying 

A e A =^ Vi (t>i{X) = I, 
h(^A=^3i (j)i{h) ^ I. (116) 

A probabilistic projection onto A can be constructed us- 
ing controlled conjugations on an ancilla \b): 

J2 ah\hbh-^) ® \b) 

heH 

ah\hbh-') (g, \(j>,{h)bMh)-^), (117) 

heH 
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and then using fusion to make sure that the ancilla re- 
mains in the |6) state. Repeating the procedure for each 
produces the desired projection. 

To build the set of functions that are balanced on A, 
let $ be the set of functions in the form of Eq. H115() 
such that (j){a) e A — {/}. We shall prove that this is the 
desired set of functions. 

Let A £ A be non-trivial and let be any map in 
$. The value of (f>{X) must be non-trivial and contained 
in A. Otherwise, it would be possible to construct a 
map in product form such that a is in its kernel but 
A is not, contrary to the definition of A. In fact, the 
functions in $ are just automorphisms of A and form a 
group with multiplication given by function composition. 
Furthermore, because (a) — H, for any non-trivial A G 
A there exists a function ipa^x € ^ such that (f>a^\{a) = 
A. If A' S A is a third non-trivial element, then for every 
function G $ such that 0(A) ~ A' there is a function 
(/)'(a) ~ A' given by 0' — cj) o <pa^\- Therefore, <I> is 
balanced on A. 



2. The amplitudes F^^^ 

To choose R and 7 we first need to examine and de- 
fine Fh^j more carefully. Since we are mostly interested 
in whether F^^-y is zero or non-zero, we will generally 
work with its magnitude squared, which has the simple 
expression 



(118) 



where Pj is the projector onto the space that will turn 
into the representation 7 after fusion. This subspace is 
just the subspace that transforms as 7 under conjuga- 
tion. It can be projected out using the orthogonality of 
characters (and matrix entries for non-abelian represen- 
tations): 



9^G 



(119) 



where 7 is the conjugate representation. Note that the 
values of the representation 7 on g G G will be denoted 
by 7^, as a reminder that it is always a power of some 
root of unity which we shall also denote by 7. 

Combining the expressions for the projector and the 
electric charge state we obtain 



geG 



dR\G\^^ 



(120) 



where dji is the dimension of representation R. In the 
second line, the magnitude squared of the matrix is given 



by |Mp = Tr (Milff), which is equivalent to the sum of 
the magnitude squared of the entries of the matrix. 

Because H is abelian, the representation R can be di- 
agonalized on H so that the diagonal entries are one- 
dimensional representations of H. These representations 
can be labeled by an index i running along the diagonal of 
the matrices R, and described by functions lj^ : H ^ C. 
With the new notation: 



\Fu 



1 



~-y 



ghg 



(121) 



where the representation R is now implicit in the defini- 
tion of the representations {LOi\. 

Finally, let S be the stabilizer of H in G, that is, the 
subgroup of G that commutes with every element of H . 
Clearly, it is a normal subgroup of G, and H C S. Fur- 
thermore, we had argued that if qi — p then Kq. G S. 
Therefore IG/S"! is not divisible by p. 

Since the function F^^^ will be zero unless we choose 
a representation such that 7'^ 
from now on, and write 



1, we shall assume this 



IF, 



dR\G\ 



1=1 



E 

geG/s 



7J [j; 



ghg 



(122) 



We are now guaranteed that 7 corresponds to powers 
of an n*^ root of unity such that p does not divide n. 
The terms in the above expression have the form 



(123) 



where the coefficients Ci are sums of n*'' roots of unity. 
By Ref . ||7j , the expression will be zero if and only if the 
n coefficients Ci are all equal. 

Using the above notation it is easy to show two prop- 
erties of the amplitudes Fh^-y- If 7^ then 



\F, 



dR\G\ 
^ 0, 



2 ^-^ 

i=l 



E 7^ 

96G/S 



(124) 



as long as p does not divide j. Note that in general 
7^ The fact that was used above is that 

^ implies that at least two coefficients of differ- 
ent powers of uj must be different. Replacing w by a power 
of itself just permutes the coefficients Ci in Eq. (|123|l . 

The second property is easier to prove in the form of 
Eq. H121() and states that given \Fh^^ \ ^ then 



\F^ 



1 



hx—^-y I 



dR\G\- 



E 



g&G 



-a qxhxq 

7^ uj^ 



dR\G\ 



d,R 
i=l 



,ghg 



geG 



(125) 
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for any x G G. The second line involves a relabeling of 
the summation variable, whereas the third line is true 
because 7 is a group homomorphism and 7"^ is just an 
overall phase. 

Together, the two properties imply that, if 
is non-zero, then so are the amplitudes |F/i'^-y| for 
any non-trivial h' = gh^g~^. Unfortunately, even after 
adding the identity element, this set is in general not a 
group. Furthermore, it remains to be shown that the 
amplitude is non-zero for at least one h. 

3. Finding a non-zero amplitude 

It is possible to indirectly show that, for every clement 
h £ H, there is a pair of representations R and 7 meeting 
our requirements, such that \Ffi^^\^ ^ 0. 

The basic idea is to consider the regular representation 
of G. Let Hq be the Hilbert space spanned by the vectors 

\9)g (126) 

for g G G. For the moment, these are just abstract vec- 
tors in a Hilbert space, and therefore we use the above 
notation to distinguish them from the anyon magnetic 
charges. 

The group G has both a left and a right action on this 
vector space, which transforms as the regular represen- 
tation in both cases. More generally, we could say that 
there is an action of the group G x G on this vector space 
given by 

I.9)g \9i992 ^)g (127) 

for any element gi x g2 G G x G. 

Let Hh be the Hilbert space spanned by the vectors of 
the form |M)^, where R in an irreducible representation 
of G. These spaces are also representations of G x G and, 
in fact, are irreducible. The space Hq decomposes as a 
sum of irreducible representations of G x G as 

^g = 0Wh (128) 

R 

with each irreducible representation R appearing exactly 
once. Fusion corresponds to a further decomposition into 
the irreducible representations of the diagonal group G. 

The state \I)q transforms as the identity under the 
diagonal group, and can therefore be written as a sum 
of states |i?(/))jj for different representations R. Hence 
a state \h)Q can be written as a sum of states \R{h))j^. 
If the state has a non-zero projection to a repre- 
sentation 7 of the diagonal group, then we know that 
|J)(^-y| ^ for at least one irreducible representation R. 

More explicitly, the projection is 

^7l^)(5 = ^E^'l5%"%- (129) 



To make it non-zero, it is sufficient to choose 7 to 
be constant over the stabilizer, Sh, in G of h. This is 
still possible, even with our requirements that 7 be one- 
dimensional and non-trivial, because Sh/H is a proper 
subgroup of the nilpotcnt group G / H. Proper subgroups 
of nilpotent groups are always contained in proper normal 
subgroups because the normalizer of the proper subgroup 
is always a larger group (and eventually the operation of 
replacing a subgroup with its normalizer must yield a 
normal subgroup). This concludes the proof that, for 
any non-trivial h G H, there exists a choice of 7 and R 
such that ^ 0. 

In fact, for any two non-trivial elements Ai , A2 G A, the 
same representation 7 is useful because Sx-^ — Sx^. How- 
ever, it is not clear that it is possible to pick R such that 
both I^Ai^tI and |-Fa2^7I ^ 0- This is illustrated by 
working with the group Zg Xg (Z2 x Z3), where certain 
choices of 7 consistent with the above discussion lead to 
zero amplitudes for at least one non- trivial element of A, 
no matter which R is used. On the other hand, the same 
example does have simultaneous choices of R and 7 that 
satisfy all our requirements. It is unclear to the author 
whether it is possible, for any group G, to choose R and 
7 such that |Ja^7| for all non- trivial elements A G A 
simultaneously. 

4- Alternative A 

What happens if R and 7 cannot be chosen so that 
|Fx^-y| ^ over all non- trivial elements A G A? While 
none of the examples in this paper have this problem, 
if such a case arises, we could try to shrink A. In par- 
ticular, if A = Zp, then the problem is solved. That is, 
because we can always choose the representations so that 
the amplitude is non-zero for some element, and then it is 
guaranteed to be non-zero for the powers of that element 
as well. 

The set of functions balanced on A = Zp can be easily 
constructed as simply (j){g) = for < z < p. However, 
the probabilistic projection onto A is more difficult. It 
can be achieved if we are willing to relax the error model 
of the probabilistic projections. That is, we use an ap- 
proximate probabilistic projection, where the probabili- 
ties and projected states are close to the desired results. 
While the results will be exponentially close in the num- 
ber of successful fusions, they will only be polynomially 
close in the number of actual fusions, and therefore the 
machinery of fault tolerant quantum computation must 
be employed. Computation with the approximate gate 
will still be feasible, but one of the advantages of topo- 
logical quantum computation, that is, the exactness of 
gates, will be lost. 

To construct this approximate projection, consider the 
amplitude for the fusion of the electric charges into the 
vacuum, denoted by Fh^i. It is the same quantity that 
has been dealt with thus far, only with the representation 
7 replaced by the identity representation. These quanti- 
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ties have the expression 
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where CQ{h) is the conjugacy class of h in G. The am- 
phtudes satisfy the properties 



< IF, 
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for any non-trivial h Cz H. The first inequality comes 
from the fact that we are summing p*^ roots of unity 
and the number of summands is not divisible by p. The 
second inequality comes from the fact that ui'^ must be 
non-constant over the conjugacy class of h. The equation 



n 



(132) 
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is true because the right-hand side commutes with all of 

G, and therefore must be the identity. Because the num- 
ber of factors on the right is not divisible by p, uji cannot 
be constant over the conjugacy class unless it is the iden- 
tity. Furthermore, since the conjugacy class generates 

H , and R is non-trivial, one of the uji must not be the 
identity. This proves the second inequality of Eq. (|131|l . 

The standard procedure of entangling a state with an 
electric charge pair, which is then fused, can then be 
used. The state is now kept if the pair fuses into the 
vacuum, which always has a non-zero probability of oc- 
curring. The basis state that was entangled with |i?(/))^ 
will have its amplitude increased relative to the other ba- 
sis states. Using braiding to achieve a function of the 
form f{h) = /la*, for some element a € A and different 
values of i, we can make the basis states in A consist- 
ing of powers of a have an arbitrarily large amplitude 
relative to the other states. Even if |F/i^/| varies sig- 
nificantly over the non-trivial elements of H , we can use 
the old A projector and functions in $ to balance out the 
non-trivial elements while increasing the amplitude of the 
state with f{h) — I. After many repetitions, the basis 
states with flux a^ba~^ can be made to have an amplitude 
much larger then all the other states. This completes the 
construction of the approximate probabilistic projection 
onto the new A for the special cases when we require 
A = Z„. 



and |0, ...,d). Furthermore, there exists a non-trivial 
subgroup A C if, such that we can implement proba- 
bihstic projections onto A and |0, . . . , O)""" fl A. 

To define the real computational subspace, choose a 
non-trivial element a g A, and define 



Wba 



(133) 



for Q < i < p. This subspace corresponds to the subgroup 
{a*} C A of powers of a. 

A probabilistic projection onto the real computational 
space, corresponding to {a*}, can be achieved in two 
steps. The first step is to apply the probabilistic projec- 
tion onto A. The second step is repeated for each A G A 
that is not in {a*}. For fixed A, we use an ancilla to con- 
jugate by A^^, then do the probabilistic projection onto 
|0, . . . , 0)^ n A and then conjugate by A using another 
ancilla: 

ax\xbx~^) — > ax\Xxbx~^ \~^) 

— > G ^ aa,\\xbx-^\-^) 
— > G ^ a^\xbx-^), (134) 

x^A,x^\ 

where the probabilistic projection was assumed to suc- 
ceed in the second step, and therefore the state is renor- 
malized by the constant G. The net effect of one such 
operation is to project out the state |A6A~^). If all the 
projections succeed, then we will have projected the orig- 
inal state into the computational basis, completing the 
probabilistic projection onto {a*}. 

For the case of qudits with d — p > 2 we are now done. 
The generalized controUed-X behaves as a controlled- AT 
when restricted to act on the computational space. A 
probabilistic projection onto |0) is just the probabilistic 
projection onto |0, . . . , 0) because |0, . . . , 0) = |0). The 
probabilistic projection onto |0, . . . , 0) behaves as a prob- 
abilistic projection onto |0) because 



(i|0, . . . ,0) cx Si^o 



(135) 



with the caveat that we must use the projection onto the 
computational basis to turn the |0, . . . , 0) ancillas into 
|0) ancillas. Finally, the probabilistic projection onto 
|0, . . . , 0)^ n A, reduces to a probabiHstic projection onto 
|0)^ when acting on states in the computational sub- 
space. These are the gates that were proven universal for 
quantum computation in Section If VI 



G. Putting it all together 

At this point we have shown the existence of an ex- 
tended computational space, with elements labeled by 
H — Zp, on which we can perform the generalized 
controlled- A", and probabilistic projections onto |0, . . . , 0) 



1. The case p = 2 

Special treatment must be given to the case when 
p = 2, that is, when working with qubits. Though all the 
gates constructed above are valid forp = 2, the gate-set is 



26 



not universal. The problem is that the probabilistic pro- 
jection onto jo)"'" = |1) does not provide any additional 
computational power beyond the probabilistic projection 

onto |0). 

Just as in Section HVBI the gate-set can be made uni- 
versal given a supply of the magic states: 

hi 

\M = ^E^'-"'^'^K)^b')' (136) 

where the second state can be produced from the first 
one by measuring the third qudit in the X basis. 

The production of the magic state |0a/i) is the step 
that requires a projection constructed from the fusion of 
electric charges. Given our choice of a g A above, assume 
that bab^^ £ A. This must be the case if A was defined 
as the intersection of kernels of functions. Clearly, we 
can apply a controlled conjugation by a, and therefore, 
additionally, the controlled conjugation by bab~'^ and by 
abab~^. Note that a ^ bab~^ because b was chosen to 
not commute with a. 

We begin with the |0) (g) |0) (g) |0) state and append a 
|0) = \b) ancilla. We then conjugate it to obtain 

111 

7^ E E E I*) ^ b') ^ 1^) ^ \k3,kbf-]^k)^ (137) 

v ° 1=0 3=0 fc=0 

where 

(138) 

with X to be determined in a moment. A probabilistic 
projection onto |0, . . . , 0)^ is then applied to the last an- 
cilla, and the conjugations are undone. 

If the projection succeeds we will have projected out 
two out of the initial eight basis states, depending on 
the value oi x £ {a,bab~^,abab~^}. In all cases, the 
state |1) g) |1) (g) |0) is removed, and for each of the three 
values of x, one of the other undesirable basis states is 
removed. Repeating the procedure once for each value of 
X produces the desired magic state \(j)Mi)- Note that the 
above procedure succeeds because = 1, and a com- 
mutes with bab~^. 

What happens if bab^^ is not in A? This is the case 
when fusions of electric charges into the vacuum must be 
used. In particular, instead of projecting out the undesir- 
able basis states, we increase the amplitude of the desired 
basis states, and obtain an ancilla that is exponentially 
close to the desired magic state. The procedure is almost 
unchanged, except that the function involved is 

f,,,^k^a'{bab-y x'-'', (139) 

and the function f^j^. = a* (bab^^Y must also be used 
to adjust the relative amplitude of |1) (g) |1) (g) |1) with 
respect to the other desired states. 



In either case, we have now shown that case of qubits 
can be dealt with in a similar fashion to the general qudit 
case, and therefore, we have completed the construction 
of universal quantum computation for anyons based on 
solvable non-nilpotent groups. 



VII. LEAKAGE CORRECTION 

Before concluding this paper, it is important to ad- 
dress the issue of fault tolerance. A physical system 
with anyons will have sources of errors due to the fi- 
nite separation of anyons and non-zero temperature (see 
Refs. 1, 3, 4] for details). While the probability of error is 
exponentially small in the distance and temperature, it is 
in general non-zero. These errors could be especially rel- 
evant if anyons are used as long term quantum memory, 
in which case error correcting codes must be employed. 

While most of the machinery of error correcting codes 
can be applied directly to anyons, it requires that states 
with errors remain within the computational subspace 
(that is, the subspace on which universal quantum com- 
putation can be done). For our model of computation, 
this is only a small subspace corresponding to anyons that 
are magnetic charges with fluxes such as a*6a^% and ar- 
ranged in pairs of trivial total flux. Note that only the 
magnetic charges need error correction as they are the 
ones in which the quantum state is stored. 

All that is required to perform quantum error correc- 
tion is to be able to replace qudits that have "leaked out" 
of the computational subspace with arbitrary states that 
are in the computational subspace. This step can then be 
followed by the standard error correcting step, which will 
remove the errors. The leakage correction step is equiva- 
lent to the swap-if-leaked gate described by Kempe et al. 

m 

In Ref. 01 a leakage correction scheme was presented 
for non-solvable anyons. While a similar scheme could 
be constructed for the solvable anyons discussed in the 
present paper, it will be easier to present a generic leakage 
correction scheme that can also be applied to anyons. 

The scheme is simply to teleport a computational qudit 
to a fresh qudit. The standard steps, shown in Figure ^ 
are first to create the entangled ancilla |$) = \i) ® 
\i) / \/d, and then measure the computational qudit and 
the first ancilla qudit in the basis |a,6) = X'^Z^ ® /|$), 
obtaining outcome a, b. The correction gate X'^Z^ is then 
applied to the second ancilla qudit, which now becomes 
part of the computational space. All these operations 
can be performed using the anyon gates discussed so far. 

If the original qudit was in the computational space, 
then its state will be flawlessly transfered into the new 
qudit (in our case, a fresh anyon pair). However, if the 
original qudit had leaked, then the new qudit will be 
guaranteed to be in the computational subspace, because 
it was obtained by applying Pauli operators to a qudit 
known to be in the computational subspace. This is the 
desired leakage correction protocol. 
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FIG. 1: Leakage correction circuit 



In fact, this scheme can be apphed to almost any sys- 
tem, as long as we can guarantee that the measurement 
of the first two qudits will not affect the third qudit in 
any way, as should be the case if they are sufficiently 
separated. 

The leakage correction scheme has caveat from a the- 
oretical standpoint, though. We are effectively assum- 
ing that we possess a classical leakage detection machine, 
through which the data "a, 5" is run. That is, if the mea- 
surement produced an outcome in the form of a voltage, 
and then the gate X°'Z^ was constructed as a Hamilto- 
nian controlled by this voltage, we would need to guar- 
antee that only the (P acceptable voltage signals could 
reach the machine operating on the third qudit. How- 
ever, in practice, leakage correcting a classical signal is 
trivial, as classical information can be measured without 
any negative side effects. 

A very similar scheme can be produced given a quan- 
tum system that is known to have exactly d states. The 
qudit is simply swapped into the new system, the first 
system is then erased and restored into the computa- 
tional space, and then the qudit is swapped back. In this 
context, the teleportation scheme is in effect a way of 
swapping a qudit into a classical system. 

Though the leakage correction scheme was discussed in 
general terms, it clearly applies to the anyons discussed 
in this paper, and its use allows quantum error correction 
and fault tolerance to be employed. We have therefore 
shown that even in the presence of small sources of noise, 
the anyons can still be used for universal quantum com- 
putation. 



The elementary unitaries are always perfect, whereas the 
measurements are either perfect, or are known to have 
failed (i.e., when none of the probabilistic projections 
succeed). This occurs with a probability that can be 
made exponentially small in the number of fusions. Of 
course, a real system will have additional exponentially 
small errors due to finite size and temperature effects. 

The physical requirements for the constructions in this 
paper include a supply of electric charge ancillas, in ad- 
dition to the requirements of Ref. The necessity of 
the electric charges may present an extra source of diffi- 
culties for a real implementation. The exception is 5*3, in 
which case only magnetic charges are required, as men- 
tioned at the end of Section ITTll In either case, the issue 
of producing the elementary electric or magnetic ancillas 
is not addressed in this paper, though a generalization of 
the construction in Ref. [l| may be sufficient. 

Another open question is whether anyons from non- 
abelian nilpotent groups are capable of universal quan- 
tum computation. Additionally, not much is known 
about computing with anyons that do not belong to 
the electric and magnetic charge model discussed in this 
paper. On the other hand, the universality of anyons 
from certain continuous groups has been discussed in 
Refs. El El. 

Of course, the most important open question is 
whether we can find a laboratory system with anyons 
out of which a quantum computer can be built. The re- 
quirement of a two-dimensional space severely limits the 
possibilities. However, certain exotic systems such as the 
fractional levels of the quantum Hall effect may contain 
non-abelian anyons. Another option is the possibility 
of engineering a system with the desired anyons. Recent 
proposals include using optical lattices E3 or Josephson- 
j unction arrays In the latter case, an explicit array is 
constructed that simulates 5*3 gauge theory on a lattice. 
Ideally, one day such a system could be used to turn the 
ideas presented here into a working quantum computer. 
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VIII. CONCLUDING REMARKS 

The main result of this paper is that anyons from fi- 
nite groups that are solvable but not nilpotent are capa- 
ble of universal quantum computation. This set includes 
many groups of small size, which are more likely to be 
found in a physical system. Combined with the results of 
Ref. T| , we have proven that every finite group that is not 
nilpotent produces anyons capable of universal quantum 
computation. 

Furthermore, except for the groups where the methods 
of Section IVI F 41 must be used, the computations with 
anyons can be made error free in the following sense: 
in the theoretical limit of zero temperature and infinite 
separation between anyons, an arbitrarily long calcula- 
tion can proceed without the need of error correction. 
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